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Multiple  scattering  of  sound  by  correlated  monolayersa> 
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Earlier  results  for  scattering  by  random  uncorrelated  planar  distributions,  and  by  doubly  periodic 
planar  configurations  of  relatively  arbitrary  obstacles,  are  generalized  to  pair-correlated 
nonsymmetrical  monolayers.  The  existing  development  for  parallel  cylinders  in  terms  of  the 
Zemike-Prins  one-dimensional  pair  function  p(x),  is  extended  to  analogous  two-dimensional 
distributions  specified  by  /?(R)  for  aligned  impenetrable  disks.  We  obtain  the  average  multiple 
scattered  transmitted  and  reflected  waves,  and  an  energy  conserving  approximation  of  the 
differential  scattering  cross  section  per  unit  area.  Simplified  forms  are  developed  to  facilitate 
determing  p  by  inverting  measured  data.  Closed  form  low-frequency  results  are  derived  for 
identical  ellipsoids  aligned  nonsymmetrically  to  the  plane  of  centers,  and  the  array  multi] 
coupling  processes  are  discussed  in  terms  of  functions  of  p  and  their  approximations.  A • 
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INTRODUCTION 

In  previous  papers1-5  we  considered  scattering  of  waves 
by  random  and  by  periodic  planar  distributions  of  relatively 
arbitrary  obstacles.  The  developments  emphasized  scatterer 
shapes  symmetrical  to  the  plane  of  centers,  e.g.,  radially 
symmetric  obstacles,  parallel  elliptic  cylinders,  or  ellipsoids 
with  a  principal  semidiameter  perpendicular  to  the  plane  of 
centers,  or  more  generally,  a  monolayer  of  obstacles  with 
reflection  symmetry  in  a  midplane  (henceforth  a  symmetri¬ 
cal  monolayer).  For  uncorrelated  random  distributions,  in 
the  course  of  development  of  a  model  for  rough  surfaces, 1  we 
obtained  the  average  multiple  scattered  transmitted  and  re¬ 
flected  waves,  and  an  energy  conserving  approximation  for 
the  differential  scattering  cross  section  per  unit  area.  For 
periodic  monolayers  (gratings3  of  equally  spaced  parallel 
cylinders  in  two  dimensions  and  the  doubly  periodic  planar 
array4  of  bounded  obstacles  in  three  dimensions),  we  derived 
the  multiple  scattered  solutions  for  arbitrary  ratios  of  spac- 
ings  (b )  to  wavelength  (A ).  The  results  for  small  b  /A  were  also 
relevant  for  random  distributions  for  negligible  incoherent 
scattering.  In  addition,  we  considered  correlated  mono- 
layers  of  parallel  cylinders5  in  terms  of  the  Zemike-Prins 
pair  distribution  function6  p{x),  with  emphasis  on  reduction 
procedures  leading  to  the  results  for  the  uncorrelated1  and 
periodic3  limits. 

The  uncorrelated  limit  corresponds  to  a  one-dimen¬ 
sional  sparse  gas  in  statistical  mechanics,  with  minimum  se¬ 
paration  (b )  of  particle  centers  small  compared  to  average 
separation  ( b );  the  periodic  limit  (b—*b )  corresponds  to  a  one¬ 
dimensional  crystal.  Using  transform  methods  we  converted 
p[x)  to  a  residue  series5  to  obtain  more  rapidly  converging 
forms  for  Z?<6  and  b~b.  For  cylinders  of  width  la  in  the 
plane  of  centers,  we  considered  2 a<b  and  arbitrary 
kb  —  2nb  /A,  in  order  to  demonstrate  that  correlation  effects 
determined  by  b  /b=pb  =  w  (the  statistical  mechanics 
packing  factor)  accounted  for  the  full  range  of  phenomena 
from  w  x 0  for  the  sparse  gas  on  to  u?-*  1  for  the  deterministic 
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Now  we  generalize  the  earlier  results  for  cylinders  to 
multiple  scattering  by  correlated  nonsymmetrical  mono- 
layers  of  bounded  obstacles,  and  emphasize  reduction  proce¬ 
dures  for  data  inversion  purposes.  The  initial  development 
for  cylinders,5  based  on  separation  of  variables  (1933),  was 
simplified  by  a  Green’s  function  procedure  (1959);  the  pre¬ 
sent  procedure  for  bounded  obstacles  (as  well  as  for  cylin¬ 
ders)  represents  a  further  simplification.  We  specialize  the 
general  representation7  for  the  solution  IP  in  terms  of  the 
multiple  scattering  amplitude  G  of  an  obstacle  in  an  arbi¬ 
trary  configuration,  with  G  related  functionally  to  the  ampli¬ 
tude  g  in  isolation,  and  then  average.  We  consider  the  aver¬ 
age  field  ( !P  >  for  a  statistically  homogeneous  ensemble  of 
configurations  of  aligned  identical  obstacles  with  one-parti¬ 
cle  statistics  specified  by  the  average  number/?  of  particles  in 
unit  area,  and  two-particle  statistics  determined  by  /?/?(R). 
For  p,  the  distribution  function  for  the  separation  R  of 
centers  of  pairs,  we  use  general  statistical  mechanics  re¬ 
sults89  for  impenetrable  aligned  identical  disks  b(R)/2,  such 
that  the  exclusion  curve  b  (the  curve  around  one  center  that 
excludes  all  others)  determines  /?(R).  In  general,  we  assume 
that  b  and  p  have  the  same  inversion  and  reflection  symme¬ 
tries  as  an  ellipse. 

Although  even  the  radially  symmetric  special  case/?(R), 
the  radial  distribution  function  for  identical  circular  disks  is 
not  known  explicitly,  the  results  derived  for  ( IP  )  and  (|1P|2) 
are  reduced  to  facilitate  determining  p(R)  by  inverting  mea¬ 
sured  data,  essentially  as  in  the  corresponding  x-ray  prob¬ 
lems.  For  the  direct  problem  of  predicting  ( IP )  for  a  particu¬ 
lar  distribution,  or  of  synthesizing  a  monolayer  with  special 
characteristics,  several  procedures  for  approximating/?  exist 
for  special  ranges  of  the  parameters.  The  leading  terms  of  the 
virial  expansion  of  p  suffice  for  sparse  concentrations,  and 
for  arbitrary  concentration  and  low  frequencies  (small  kb )  a 
required  integral  of p  is  available  as  a  simple  function  of  the 
packing  density  in  the  plane.10  More  generally,/?!/? )  can  be 
obtained  numerically  from  the  Percus-Yevick  or  other  ap¬ 
proximations8,9  for  the  pair  function.  For  identical  aligned 
ellipsoidal  particles  that  intersect  the  plane  of  centers  as  an 
ellipse  a(R),  the  exclusion  curve  is  the  symmetrical  ellipse 
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2a(R).  For  greater  generality  (essentially  as  for  the  corre¬ 
sponding  volume  distributions10)  we  take  the  exclusion 
curve  b{R)  as  neither  similar  nor  similarly  aligned  with  a(R). 

In  the  following,  for  brevity,  we  use  ( 1 :9)  for  Eq.  (9)  of 
Ref.  1,  etc.  Section  I  introduces  required  notation  and  repre¬ 
sentations  in  the  form  of  a  brief  sketch  of  the  single  and  many 
body  scattering  problems,  and  Sec.  II  derives  ( V )  for  an 
incident  plane  wave  <f>,  and  the  functional  equation  for  the 
average  multiple  scattered  amplitude  G  [g].  Section  III  con¬ 
siders  the  average  energy  functions,  and  provides  physical 
motivation  for  Sec.  IV  in  which  G  is  expressed  in  terms  of  a 
transform  of  g  that  delineates  the  loss  mechanisms  in  the 
array,  and  that  leads  directly  to  simplified  results  for  data 
inversion.  For  the  nonsymmetrical  monolayer,  for  incident 
waves  <f>  and  <f> '  that  are  images  in  the  plane  of  centers,  we 
develop  explicit  approximations  in  terms  of  p  and  g  for  the 
magnitudes  and  phases  of  the  transmission  and  reflection 
coefficients  and  for  the  differential  scattering  cross  section 
per  unit  area.  Section  V  uses  Fourier  series  in  G  [g]  to  obtain 
algebraic  systems  for  the  corresponding  multiple  scattered 
multipole  coefficients  A  [a]  in  terms  of  the  distribution  inte¬ 
grals  that  characterize  the  array  (i.e.,  in  terms  of  contin- 
uium  analogs  of  the  lattice  sums  for  the  periodic  cases).3,4  We 
consider  3P  in  terms  of  p  (virial  expansions  and  Fourier 
transforms),  and  develop  low-frequency  approximations. 
Section  VI  derives  closed  form  results  for  small  ellipsoids 
aligned  nonsymmetrically  to  the  plane  of  centers,  with  em¬ 
phasis  on  multipole  coupling  effects.  These  explicit  results 
provide  illustrations  of  the  general  considerations  of  Sec.  IV 
for  <f>  and  <p '  incident  on  nonsymmetrical  monolayers.  For 
small  kb,  except  for  sparse  packings,  wavelength-indepen¬ 
dent  multipole  coupling  effects  may  be  misinterpreted  as 
changes  in  the  shape  of  the  obstacles;  e.g.,  because  of  array 
coupling,  a  spherical  obstacle  corresponds  to  an  equivalent 
ellipsoid  flattened  along  the  array  normal  and  broadened  in 
the  plane  of  centers. 

I.  REPRESENTATIONS 

Consider  a  plane  wave  <pe  ~  io“  incident  on  an  obstacle  (in 
two  or  three  dimensions)  with  center  at  r  =  0,  the  center  of 
its  smallest  circumscribing  circle  or  sphere  (of  radius  a).  We 
take 

*  =  «*V=/f,  r*  =  z2  +  R2 

?  =  )  =  i  cos  0  +  R(g> )  sin  0, 

A 

R(^> )  =  St  cos  q>  -I-  y  sin  <p; 

k  =  kk,  k  =  2v/A, 

k  =  H0o,<po)  =  f0  =  2  cos  0O  +  Rq  sin  0„, 

Ro  =  R(*>o)-  (1) 

When  convenient,  we  work  with  direction  cosines 

f  =  ?7  +  £a +  =  r  =  Y\  =  cos0, 

a  =  y2  =  sin  6  cos  <p. 

P  =  Yi  =  sin  6  sin  <p. 

For  cylinders  with  axes  along  y,  we  set  <p  =  0,  P  —  y3  =  0.  In 
general,  we  use  terminology  and  procedures  that  apply  for 
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three  dimensions  and  can  be  specialized  to  two,  and  treat 
two-  and  three-dimensional  problems  in  parallel;  if  two 
forms  of  a  function  or  operator  arise,  we  list  the  cylindrical 
first. 

In  the  region  external  to  the  scatterer  (outside  the  vol¬ 
ume  y  bounded  by  >¥),  the  field  ip  =  (f>  +  u  satisfies 
(V3  +  k  2)ip  =  0  with  u  as  a  radiative  function7 

«(«•)  =  c0J [k0(k  |r  -  r'\)dnu[r')  -  udnh0]dy(r')={h0,u\- 

(2) 

h0{x)  =  Hfo\x),  h (0n(x);  c0  =  1//4,  k  /iAv. 

Here,  y  is  the  obstacle’s  surface  (or  any  surface  that  isolates 
it  from  the  field  point  r),  and  d„  is  the  outward  normal  deri¬ 
vative.  For  r~  oo , 

u~h  (kr)g(f,k),  g(#,k)=  {e‘k'  r',a},  k,  =  kr,  (3) 

with  h  as  the  asymptotic  form  of  h „,  and  g(f,k)  as  the  scatter¬ 
ing  amplitude.  Substituting  the  plane-wave  spectral  repre¬ 
sentation  for  h0  in  (2)  we  obtain,  at  least  for  r  >  a  for  all  ?, 

a(f)  =  Je*'rg(fc,k), 

kc=krc,  ic=r{0c,<pc),  (4) 

where  Sc  equals  (l/ir)J  d&c  with  contour  as  fori/),",  or 

with  contours  as  for  h  j}1.  Decomposing  g  as  Fourier  series 
(trigonometric  or  spherical  harmonic)  of  scattering  coeffi¬ 
cients  a„(?o)  or  a”(?0),  we  obtain7  from  (4)  for  u  the  corre¬ 
sponding  series  in  special  functions  H[l\kr)e>"e  or 
At„,,(kr)T”(f).  We  reserve  such  series  for  special  computa¬ 
tions,  and  use  primarily  (2)  and  (4)  to  derive  general  results. 

The  integral  in  (4)  corresponds  to  h0.  In  terms  of  the 
associated  integral  operator  for  j0  =  Re  the  mean  (JY) 
over  all  real  directions  of  observation,  we  write  the  energy 
theorem  as 

—  Go  Reg(k,k)  =  a„  +as; 

°o  =  -\a0{rp*,ip\,  a,  =  \a0[u*,u]  =  tv*|g<f,k)|2;. 

(5) 

_  *  4tt. 

<Tf>~T'~kji 

=  frdd,  —  Cd<p  r dd  sin  0, 

2irJo  4 it  Jo  Jo 

where  oa  and  a,  are  the  absorption  and  scattering  cross  sec¬ 
tions.  This  theorem  follows  directly  from  the  definition  of  <xa 
in  terms  oiip  =  <P  +  u  by  using  Green’s  theorem  and  (3).  The 
same  procedure  applied  to  $?,)  and  $?2),  the  solutions  for 
two  different  directions  of  incidence  t  „f2  (each  solution  sa¬ 
tisfying  the  same  conditions  on  y  and  in  V),  gives  the  usual 
reciprocity  relation 

{ =  p:  gi  -  f  „?2)  =  g<  -  f2,f .), 
g(f,k)=g(-k,t-f).  (6) 
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The  following  development  applies  for  general  obstacles  and 
all  conditions  on  y  and  in  for  which  (6)  holds. 

Fora  fixed  planar  configuration  of  N  obstacles  with 
centers  (R, )  in  the  xy  plane,  we  write  the  solution  external  to 
all  obstacles  as 

*  =  <f>  +  f  U,(t  -  Rj;R„R2)...,R„ )***'  =  </>+&, 

5=  1 

R,  =xJi+^y  =  RJR(^J), 

C/5(r-RJ)  =  J^r-R*'GJ(fc,k), 

Gs(r,k)  =  {e-^',Us\„  (7) 

where  the  multiple  scattering  amplitude  Gs  is  the  form  (3) 
over  ys(t'),  and  r'  is  the  local  vector  from  R,.  For 
|r  —  R,  |  ~  oo ,  we  have  Us  ~h  ( k  |r  —  Rs  |)GS.  With  reference 
to  the  center  of  scatterer  t,  we  display  the  phase  factor  intro¬ 
duced  by  </>  (R, )  and  use 

0  =  V,e*R'  =  (0,  + 

0,  =  e^'  +  2; US(H,  +  r'  -  R>*‘,R'  -  "■>, 

where  2'  is  the  sum  over  s^t.  We  require  that  .0,,  0,,U, 
satisfy  the  same  conditions  on  y,  and  V ,  as  for  the 
corresponding  scatterer  in  isolation,  and  that 
(V2  +  k2)U,  =  0  external  to  2^,.  Consequently  from 
{ 0(?,),0(r2)J,  =  0,  we  have 

G,{  -  [0,(?2), «(?,)),; 

substituting  for  0,  with  U,  as  in  (7),  and  using  the  definitions 
of  g  and  G,  it  follows  that 

<?,(?,£)  =  g,(?,k)  +  R-g,(f,fc)<?,(fc,k),Rtt 

=  r,  ~  rJ>  (8) 

at  least  if  the  scatterers’s  projections  do  not  overlap  on 
R„  =  R„/R„.  The  functional  equation  (8)  is  essentially  a 
reciprocity  relation  between  the  multiple  (G )  and  single  (g) 
scattering  amplitudes,  for  both  the  direct  and  inverse  scat¬ 
tering  problems.  We  emphasize  the  direct  problem  G  [g]  as  if 
g  were  known,  but  the  results  also  apply  for  the  inverse  g[G  ] . 

II.  THE  AVERAGE  WAVE 

We  average  0  over  a  statistically  homogeneous  ensem¬ 
ble  of  planar  configurations  of  identical  aligned  obstacles 
with  one-particle  statistics  specified  by  the  average  number 
ip)  of  centers  per  unit  area.  Thus 

<0 )  =  ^  +pjelk'n-  <Us(r  —  R,)),  =<f>  +  {<%),  (9) 

where  (U,),,  the  average  with  one  variable  held  fixed,  de¬ 
pends  only  on  R,  (now  a  dummy).  From  (7), 

<£/,(? -R,)>,  =J>(r-R-><G,(fc,k)>, 

=  JyM"R-,G(fc,k),  (10) 

with  (G,),  =  G  independent  of  R,,  because  was  fac¬ 
tored  at  the  start. 


Substituting  (10)  into  (9),  we  write  R,  =  +ijy, 
d R,  =  didr,,  and  (f0  -  fc)-R  =  (a0  -  +  ifio-0c)v, 

the  integral  over  £  reduces  to 

J"  ft'**.-**.  2trg(a0-ac) 

and  that  over  r,  to  2ir8{J30  —  0C  )/k.  The  8  functions  cor  re  - 
spond  to  ?c  =  r(G0,<p0),r(7j-  —  90,<p0)  =  =  k,k'  forz^O. 

Transforming  $d{ic  to  SSdac  dfic  y~  1  in  terms  of  real  varia¬ 
bles  acfic  [with  ye  =  (1  -  a2  -02)I/2  =  \yc\  or  i\yc\  for 
a2  +  P 2  <  1  or  >  1],  we  reduce  the  average  scattered  wave 
(  ^ )  as  before1  to  the  specular  values 

<^>>=*2CG(k,k)^2CG00, 

{*%  <)  =4>  '2CG  (k',k)  =  <f>  '2CGV0; 

C  =  p/ky0t  prr/k2y0;  <f> '  =  e*-', 

k'  =  *o  -  ?(«■  -  &o><Po)  =  k  -  2k££,  (11) 

where  <f> '  is  the  image  of<^  in  the  plane  z  =  0.  The  subscripts 
( >  and  < )  indicate  that  the  forms  hold  at  least  for  z  >  a  and 
z  <  —  a,  respectively.  (Were  the  incident  wave  <f> '  instead  of  <f> 
then  k  and  k'  would  be  interchanged.)  The  corresponding  net 
average  fields  in  the  transmitted  and  reflected  regions  are 

<0>>=^(l  +  2  CG00)  =  <f>T, 

<0<>=^  +  ^'2CG„o  =*  +  f'R,  (12) 

with  T  and  R  as  the  coherent  transmission  and  reflection 
coefficients. 

The  two  values  in  (11)  correspond  to  one  wave  form 
multiplied  by  the  forward  scattered  (G,*,)  or  reflected  (G„.0) 
amplitudes: 

(%)  =e± ikzr° +  **2CG (  ±  2  cos  0O  +  Ro sin  ^ko) 
a  e‘k  \z\r«  + 

is  a  wave  outgoing  from  z  —  0.  The  images  <t>  and  <f> '  have  the 
same  \z\pcj>  dependence  and  comprise  the  central  propagat¬ 
ing  mode  of  the  periodic  cases3-4  discussed  earlier;  the  earlier 
results  for  the  case  of  one  propagating  mode, 
kb  ( 1  ±  sin  90)  <  2 it,  provide  bounds  for  present  low-frequen¬ 
cy  approximations. 

Similarly,  <G(>,  =  G  is  obtained  from  the  ensemble 
average  of  (8)  for  pair-correlated  scatterers  with  separation 
(Ru)  specified  by  pj>(„ ),  such  that  p( R)~  1  for  R~oo,  and 
p(R)  =  Ofor/?  <  b  (R)  with  A  (R)as  the  minimum  separation  of 
pairs.  We  assume  that  p(R)  is  fully  determined  by  the  curve 
R  =  b(R),  the  exclusion  curve  enclosing  the  area  containing 
the  centerof  only  one  scatterer.  For  radially  symmetric  sta¬ 
tistics,  b  (R)  =  b  is  constant  and  p  =  p{R  )  is  independent  of 
direction.  More  generally,  we  consider  exclusion  areas  hav¬ 
ing  the  inversion  and  reflection  symmetries  of  an  ellipse, 
henceforth  elliptical  symmetry.  From  (8), 

<G, (?,£)),  =  g(r,k)  +p  JdR,  p(Ru)e -*«'■ 

xj*(Mc)«7,(  fc,k)>„AR“,  (13) 

where  (G, )  „  is  the  average  over  all  variables  but  Rs  and  R, . 
To  obtain  a  deterministic  equation  for  (G, ),  =  G,  we  pro- 
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ceed  as  before1,5  and  replace  {G,)s,  by  (Gs)s  =  G.  See  dis¬ 
cussion  after  (10:16)  for  citations  to  work  of  Lax  on  the  corre¬ 
sponding  approximation  for  random  media,  and  for 
alternatives  by  Keller  and  by  the  writer. 

We  consider 


+pjdRp(R)e  '  «J>*  -  V-fc)G(fc,k); 

(14 

HA 

I  /Mr/2  —  i  OO  |  pv/2  —  loo 

=  -|  d$c,  dip  A  dQc  sin  6C , 

IT  J  —  ir/2  +  i  oo  2lT  J  —  tt  Jo 


in  terms  of  simple  contours,  and  z  =  ze  to  provide  a  conver¬ 
gence  factor  with  exponent  ikz  cos  6C  =k  |e|sinh  r  for 
9C  =  ±  (ir/2  —  it)  and  r  ~  oo .  We  asume  that  G  (?,k)  may  be 
equated  with  the  limit  of  (14)  for  z=|c|— >0,  or  for 
z—  —  |e|— >0  and  fc  replaced  by  its  image  i ' ;  in  general,  we 
work  with  the  mean  of  the  two  which  makes  the  symmetry  of 
the  problem  more  evident.  The  form  &r(tc)  =  &(tc) 
+  #(?'),  with  #(rc)  =  g(r,rJG(?<.k),  depends  only  on  the 
real  variables  <pc  and  sin0c  (i.e.,  sin  9C  =  ^  cosh  r  is  real 
even  for  9C  =  ±  rr/2  :f  ir)  as  may  be  seen  by  using  trigono¬ 
metric  or  spherical  harmonic  series  representations.  For  cyl¬ 
inders,  ^(t)  depends  on  exp[//i  0]  +  exp[/n(w  —  6 )],  so  that 
only  cos  2nd  and  sin  (In  +  1)0  arise;  these  functions  are  pol¬ 
ynomials  in  sin  0.  Similarly  for  bounded  obstacles, 

F;(cos^H-/,"(-cos0)=  [l  +  (-  l)',  +  'n]f,:r(cos0) 

selects  Pl7  and  and  these  too  are  polynomials  in 

sin  0.  The  exponent  k  (fc  —  f0)-R  =  k  (?'  —  r0)R=KcR 
also  depends  only  on  sin  6C  and  <pc,  so  that  ekc'RJr(r<.)  is 
symmetrical  to  reflection  in  z  =  0. 

We  introduce  the  Fourier  transform  of  the  pair  correla¬ 
tion  function 


P[KC]  =p\dYlp(R)eX'*, 

(15) 

Kf  =-  (ke  -  ko)-(ti  +  ff)  =  k  (ac  -  a0)i  +  k(/3c-  P0)y, 


a  real  function  of  Kc  for  elliptical  symmetry,  and  briefer 
notation  G (f,?0)  =  GM,G(fc,i0)  =  G^,  etc.  Thus  the  mean 
form  of  (14)  may  be  written  essentially  as  (5:36). 

G^-g*,  +  +gn’GC’0) 

S  =  lim—  [e~k^P[KA, 

*—0  2  Je 

S  =  S,  +  S„  Sr  =  Re S  =  ,P [Ke].  (16) 

As  before  for  the  periodic  cases,  S,  corresponds  to  real  0C 
(the  propagating  range),  and  S€  to  \9C  \  >  ir/2  (the  evanescent 
range).  The  operator  Jt ,  representing  the  mean  over  the 
right  half-space  (\9 1 <w/2)  can  be  replaced  by  of  (5)  be¬ 
cause  of  the  symmetry  of  the  operand.  The  operator  S,  is 
imaginary;  using  9C  =  ±  (w/2  —  ir)  yields  real  integrals 
times/. 

Equivalently,  in  terms  of  the  Fourier  transform  of  the 


total  correlation  function  p(R)  —  1, 

£P[K  ]  =pj [p(R)  -  1  ]e*"dR,  (17) 

which  differs  from  P[K]  by  the  5-function  contributions  at 
?0,?o » we  make  specular  terms  of  S  explicit.  Thus 

6,0  —  i*>  +  C  [gMGM  +  gtfGff o  ] 

+  Sm(gnGc0+g^G,0),  (18) 

where  St0)  differs  from  S  of  ( 1 6)  in  containing  @  [K  ]  instead 
of  P  [K  ].  If  the  effects  of  S,0)  are  negligible,  then  ( 1 8)  reduces 
to  the  form  analyzed  before  for  sparse  uncorrelated  distribu¬ 
tions.  The  earlier  construction  consisted  essentially  of  taking 
the  multiple-scattered  excitation  of  one  obstacle  in  the  array 
as  the  mean  of  the  average  transmitted  and  reflected  fields 
H<^>>  +  <^<>]°f(12). 

In  Ref.  5  for  correlated  cylinders,  we  emphasized  re¬ 
duction  procedures  of  (16)  that  led  to  earlier  results  for  the 
periodic  and  uncorrelated  cases,  but  now  we  exclude  period¬ 
ic  distributions  (except  for  the  low-frequency  case  of  one 
propagating  mode)  and  emphasize  the  effects  of  correlations. 
To  motivate  the  reduction  for  G,  we  first  consider  the  aver¬ 
ages  of  the  quadratic  function  of  the  field. 

ill.  AVERAGE  ENERGY  FUNCTIONS 

For  a  fixed  configuration,  the  normalized  average  ener¬ 
gy  density  is  given  by  |  IF  |2  and  the  corresponding  flux  den¬ 
sity  by  Re(  !F  *  V!F  /ik )  with  !Fas  in  (7).  The  ensemble  average 
of  |!F|2  equals 

(|!F|2)  =  |<!P>|2  +  V, 

V  =pf(\U,\ 2>,  dRs  +p2ff[(UsUr)„p( Ra) 

-{Us)s(UT),]e~lkn°dR1dR„  (19) 

with  |  (!F  >  |2  =  |^  +  (^  )  |2  in  terms  of  (12)  as  the  coherent 
intensity,  and  the  variance  V=  ( | IF  —  < IF >  |2>  as  the  inco¬ 
herent  or  fluctuation  intensity.  Substituting  the  spectral 
form  of  Us,  and  integrating  over  R,  to  obtain  5-functions 
which  eliminate  one  of  the  spectral  operators,  we  obtain  for 
z=\z\, 

V „  =|2Cce-2Im^(<|G,(fc,k)|2>J 

+  p  fdRe*c‘*[p(R)<G,G?>M  —  (G.XG*),]^  , 

(20) 

where  Cc  —  g/y*.  For  z~  oo,  we  have  V> 

~2CygJ(r  sec  9  (  j.  Replacing  < G,G*)„  by  <G,>,  <G*>, 
with  (GS),  =  G, 

V ,  ~2CyvMr  ]  |G (f,k)|2  sec  9  =p  f  sec  9-, 

Y  > 

CK.=c12T|C?K,|2;c1=-^,-ir;  (21) 

%T[K]  =  1  +  &(K)  =  1  +pj(p-  l^dR, 
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where  9^[K  ]  is  the  structure  factor,  and  pQM  is  the  multiple 
scattered  differential  cross  section  per  unit  area.  In  (21),  we 
integrate  f,  =  S>dd,dO  over  the  forward  half-space 
\&  |  <ir/2;  for  the  back  space,  we  replace  ?  by  or  write/.,  to 
integrate  over  the  appropriate  values  of  9. 

For  the  corresponding  average  flux,  we  write 

J  =  Re<!P*VlP/i* )  =  ¥  +  I,  ¥  =  Re<!P*>V<!P>//*, 

I*  ~pj  Q^r  sec  9\  (22) 

j;  =  I T  |2k,  =  k  +  Re(e  ~  akr<*R  )(k  +  k')  +  \R  |2k'. 

The  components  normal  to  the  distribution  satisfy 

*-Jc>  =  \T\2y0,  i¥,  =  (1  +  |R  | Vo, 

*<i>  +ij=p[J  +  J]&o=p^2H<7,oI2=p^. 

with  py,  as  the  average  incoherent  scattering  cross  section 
per  unit  area.  The  corresponding  average  absorption  cross 
section  per  unit  area  is  the  average  inward  flux  —  z-(J> 
+  J  <  )  =  pya ,  and  conservation  of  energy  is  shown  by1,3 

i  =  |r I2  +  \  t\2 +py sec Oq,  y  =  ya+y,x 

R  =  2CG„0,  T=  1  +  2CGoo-  (23) 

This  states  that  the  incident  flux  density  equals  the  flux  co¬ 
herently  reflected  and  transmitted,  and  incoherently  scat¬ 
tered  and  absorbed  by  the  area  of  distribution  irradiated  by 
unit  area  of  incident  wave. 

To  obtain  y  „  as  a  well-defined  surface  integral  of  a 
quadratic  function  of  'P,  as  well  as  additional  relations  for  G, 
we  write  the  energy  absorbed  by  a  fixed  set  of  scatterers  with¬ 
in  a  closed  surface  A  as 

-  J  Re(^P)^A  =  -  |<r0{  }„ . 

Using  Green’s  theorem  and  (V2  +  k  2)!P  =  0  outside  all  ob¬ 
stacles  we  obtain  { 'P*,'P )  =  2,  ( 1P*,<P  j,  in  order  to  replace 
A  by  the  set  of  scatterers’  surfaces  (or  any  set  of  surfaces  that 
isolates  each  scatterer  from  the  others).  Working  initially 
with  *P  =  <f>  +  ^ ,  and  then  decomposing  ab ,  we  use  Green’s 
theorem  and  (V2  +  k  2)U,  =  0  outside  scatterer  s  to  obtain 

—  2  Re{^  •,£/,),  +  ( Uf,U, }, 

+  2  Re{Uf,2'C/,e_*'R"|,. 

Multiplying  through  by  —  £  a0,  we  relate  the  absorption 
cross  section  of  obstacle  t  to  the  energy  derived  from  <p  by 
interference  with  U,,  and  to  the  scattering  cross  section 
\  a0\  ah*,ab\,  of  obstacle  t.  The  decomposition  of  ( aP*,ad\, 
in  terms  of  U,  and  U, ,  corresponds  to  an  intrinsic  compo¬ 
nent  ( U*,U,},  find  to  a  set  of  two-particle  terms.  Proceeding 
essentially  as  for  (8),  we  obtain 

f  **,*),  =  2  Re  G,(k,k)  -(-  2^|G,(f,k)|2  +  S, 

y  =  2Re2'e~  Je^-G,  (fc  ,k)G  *(f?,k),  (24) 


where  the  term  in  y  follows  from  the  asymptotic  form 
U,  ~hG,  or  from  the  Jc  form  in  (7).  For  lossless  particles 
=0,Eq.  (24)  is  a  special  case  of  (7:A  17). 

We  average  (24)  with  scatterer  t  fixed  for  a  planar  distri¬ 
bution  R(xj>)  by  proceeding  as  for  (13).  Thus,  initially  we 
write 

-  Re<G,(k,k)>,  =  y a / Cq  +  ^<|G,(f,k)|2),  +  S'; 
y  =  ReS[^„(fc)  +  3  u$'c)), 
^1,(fc)=<G1(?c,k)Gf(f?,k))„,  (25) 

where 

is  the  average  absorption  cross  section  of  obstacle  t,  and 
where  S  (in  which  the  implicit  />[AT]  is  real)  operates  on  a 
function  of  the  real  variables  sin  6C  and  <pc.  Replacing 
(G,G*)SI  by  |G  |2,  etc.,  and  using  briefer  notation  we  have 

,/-'  =  ReS[|Ge0|2+|Gc.0|2] 

=  S,  [  IG^I2  +  |G/0 12]  =yP[K]  |G,<,|2, 
and  make  the  specular  contributions  explicit  by 
y  =  [K]  IG*  |2  +  C  [  | Goo  I2  +  \G<ro  |2]. 

Thus  introducing  y  into  (25)  and  combining  the  terms 
by  using  Tf  —  1  +  & ,  we  obtain 

-  Re  Goo  =  C  [  |Goo|2  +  | Gm  |2]  +  +  y.V*9 

y,=a^W\K)\G^\\  (26) 

which  also  foUows  from  (23)  on  writing/!  and  Tin  terms  of  G 
and  using  p  sec  0q/4C  =  l/a0. 

Similarly,  corresponding  essentially  to  (7:A16),  we  have 

{^•(k,),!^)},  =  Gf(k2,k.)  +  Gf(k„k2) 

+  2^G,(f,k2)G  *(f,k,) + y  21  +  Sti , 

jr21  =  2;e  "  *,R“|AR“GJ(fc,k2)G*(f?,k1).  (27) 

For  planar  distributions  R(i,$),  and  k„k2  each  correspond¬ 
ing  to  either  k  or  k',  we  have  k,-R  =  k2-R,  and  can  proceed 
essentially  as  for  (24)  in  terms  of  the  same  P[K  ]  Thus 

-  <G,(k„k2)>,  -  <G?(k*ki)>, 

=  2^2,/<r0  +  2ur<G,(f,kJ)G*(f,k1)>,  +  y2l  +  ytf; 
Si,  -S[9l%)  +  9»M]msyv, 
SfM=(GjeckI)G*{fc,k,)>„:  yi2  =  {yi'f,  (28) 

where 

—  2y\'  =  <t0\  (!P*(k,),!P(k2)>, ),. 

Dropping  the  dependence  on  s,  t  as  before,  we  write 
9  c  =  Gc2G*,  and  ?"  =  9 e  +  9,,  to  obtain 

y2l  +y*  =  sy  +  (Sir*]*  =  (S-t-S*)Jr  =  2SrJr 

=  2~#P[K  ]Gr2G*[ 

=  2y0>[K)GnG*x 

+  2C  [  G02G  *i  +  GV1 G  *  j  ] , 
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which  reduces  (28)  to 

—  Gxi  —  G*\  =  2C  [  G02G oi  +  G0'2  G  J,  ] 

+  2^3T[K]Gr2G*l+  2yil/a0. 

(29) 

Equation  (29)  represents  four  relations,  with  (26)  as  the  case 
?,  =  ?2  =  ?0  =  h. 

Similarly  from  { <^(f,),!P(f2))t ),  =0,  we  obtain  the 
reciprocity  relation 

G  (f  „f2)  =  G  (  —  r2  —  r,),  Gt2  =  G_2_,, 

for  r ,  and  r2  each  corresponding  to  either  r0  or  . 

From  (29),  in  terms  of 

T=  1  +  2CGoo,  T '  ■-=  1  +  2CG0  0. , 

R  =  2  CGq-o,  R  =  2CG00,-> 

S,2  =  sec  0O; 

S  =  5ow  =  -So-o'  >  5  =  ^o-o  =  Sm , 
we  obtain  four  equalities 

i  =  i*i2+i7T+s=i*'i2+in2+s', 

0  =  (AT*  +  H*r  +  S,')  =  (  )♦.  (30) 

The  first  two  correspond  to  conservation  of  energy  for  inci¬ 
dent  <j>  or  (f> and  the  second  two  provide  additional  con¬ 
straints.  For  lossless  particles  and  negligible  fluctuation  ef¬ 
fects  (or  for  the  one-propagating  mode  periodic  cases),  we 
drop  the  S  terms  to  obtain  the  earlier  results  discussed  in 
detail  for  periodic 12  structures;  for  such  cases,  |T|  =  |T'|and 
|J?  |  =  |J?  '|,  and  the  sum  of  the  transmitted  phases  differ  by 
180“  from  the  sum  of  the  reflected  phases.  For  a  symmetrical 
monolayer  [R  =  R',  T=  T'),  (30)  reduces  to  two  relations 

l  =  |/J|2  +  |7’|2+5,  2ReR*T+S'=0; 
l  =  |r±R|2  +  (5±5'),  (31) 

with  S and  5 "  real.  If S is  negligible,  then  Re  R  *  T  =  0,  and  R 
and  T  are  90°  out  of  phase;  for  such  cases  \T±R  \  =  l. 

IV.  SCATTERING  AMPLITUDE  REDUCTIONS 

The  procedure  we  followed  to  obtain  the  generalized 
energy  relation  (29)  corresponding  to  G  of  (18),  may  be  ap¬ 
plied  to  a  modified  scattering  problem  shorn  of  specular 
losses  specified  by 

9ro  ~  8rO  +  S10'  [grcPtO  +  grc  9c0  ]  •  (32) 

For  lossless  scatterers,  the  procedure  yields 
9,2+9?,  =  -  2^W[K)9r29?„ 

7T[K)  =  1  +  0>{K}  =  1  +pj  [p(R)  —  1  ]e'*'H  dR, 

(33) 

where K*R=  k  { f  —  ?,)*R  =  k  (?  —  ?2)*R,  and  r ,  orr2  each  rep¬ 
resents  ■either  k  or  k'.  In  terms  of  9,  we  reduce  (18)  to 

• 

G,o  =  9, o(l  4-  CGqq)  4-  9«CG0.0,  (34) 

the  form  considered  earlier1  in  terms  ofg.  The  interpretation 
of  G  [p]  is  the  same  as  before  in  that  G  corresponds  to  9 
excited  by  the  mean  of  the  averge  transmitted  and  reflected 
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fields,  of  (12),  £  [  (  <P>  )  4-  ( !P<  )  ] ,  evaluated  at  z  =  0.  The 
following  indicates  the  physics  implicit  in  9  to  provide  a 
guide  for  developing  practical  forms  for  data  inversion. 

The  analog  of  (33)  for  a  lossless  isolated  obstacle  is 

8a  +  g?,  =  -  2 ^gngr,  =  -  2J(gu£r,  (35) 

whereg12  =  g(r  „r2)  with  r,  and  f2  arbitrary;  the  second  equa¬ 
lity  follows  from  reciprocity  and  the  fact  that 

f)  =  —  r).  The  special  value  Re  =  —  JP\g^  |2 

corresponds  to  (5)  for  oa  =  0.  To  develop  approximations  of 
g  as  well  as  of  the  multiple  scattered  amplitude  for  periodic 
distributions,3,4  we  expressed  g  before  in  terms  of  the  ampli¬ 
tude  g'  for  the  radiationless  problem 

g,2  =  g'l,  +  ^g\rgr2  =  g'n  +  ^gtrg’r 2  (36) 

such  that  for  the  lossless  case,  g |2  4-  g'2f  =  0.  The  leading 
term  of  ( 36)  includes  phase  and  absorption  effects  (that  domi¬ 
nate  for  various  special  cases)  but  no  radiative  losses. 

Decomposing  the  operator  S,D|  in  (32)  as 

+  se  =  j  -  j  jpr  +  s„ 

and  representing  the  operand  by  [gp],  we  use 
v^/(f)  =  r  [/(?)  -f  /(?')]  for  (36)  to  construct 

9*>  =g<o  +  &  +S,][g^] 

=  g^  +  [y^r^+S,][g'9). 

Suppressing  Sf  by  introducing  the  radiationless  distribution 
amplitude 

^=^+S,[g'g'],  (37) 

such  that  for  the  lossless  case  g"2  +  g2I*  =  0  with  f ,  or  r2 
each  either  k  or  k',  we  obtain 

9t2=g?2+^»'lK][g7rpr2]  (38) 

as  the  analog  of  the  isolated  obstacle  result  (36). 

For  large  kb,  the  operation  over  Sc  is  negligible  for  the 
present  problem  (which  excludes  complex  00,  as  well  as  mul¬ 
timode  periodic  cases).  We  may  also  neglect  S,  for  small  kb, 
provided  that  the  scatterer’s  width  is  small  compared  to  min¬ 
imum  separation  (a/b  small);  if  not,  then  as  shown  subse¬ 
quently,  S,  may  introduce  ^-independent  factors  arising 
from  multipole  coupling.  Neglecting  Se  corresponds  to 
9u~g,2  +^^g\r9r2  which  differs  from  gl2  in  that  the 
radiation  integral  is  modulated  by  the  structure  factor 
W  =  1  +  SP ;  the  first  term  of  W  corresponds  to  radiation 
over  the  continuum  as  for  the  isolated  case,  the  SP  term 
corresponds  to  the  effects  of  correlations,  and  their  sum  rep¬ 
resents  fluctuation  scattering.  The  similarities  in  the  forms 
for  p  and  g  indicate  how  we  may  carry  over  existing  results 
by  inspection  of  g.  In  particular,  for  various  practical  prob¬ 
lems  for  small  scatterers,  we  may  construct  approximations 
for  9  by  using  phase  and  absorption  terms  ofg  and  by  modi¬ 
fying  the  radiation  term  of  g  by  incorporation  W. 

We  reduce  (34)  by  eliminating  G  from  the  right-hand 
side,  and  write 

Gr  1  =  [?rl(!  -  G922)  +  C9r292,  ]/A 

D  —  (  1  —  Cpn)(l  —  Cp22)  —  C29l2p2, 

=1  —  C{pu  +  922)  +  G2||^||; 

Gu=(9„-C\\p\\)/D,  G2!=  9^/D,  (39) 
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where  1,2  correspond  to  0,0'  for  incident  <f>,  and  to  0',0  for 
incident  <t> '.  If  we  construct  Jt  7T Grl  G  *  (for  r ,  and  r2  differ¬ 
ent  or  equal)  and  apply  (33)  to  eliminate  integrals  of  products 
of  p  s,  we  rewrite  the  results  collectively  as 

-  Gn-Gl 

=  2 C  [G02G&  +  G0.2G$t  ]  +  2JnTGr2G*x ,  (40) 

i.e.,  as  (29)  for  y  a  =  0.  If  r,  =  f2,  we  obtain  the  lossless 
version  of  (26) 

-  Re  G„  =  C(|<?0,f2  +  | Gvl \2)  +  S%T\Gri  |2,  (41) 

where  1  equals  0  or  O'.  In  the  isolated  obstacle  result 

—  Reg,,  =M  \fcrl  |2,  the  term  —  Reg,,  corresponds  to  the 
energy  lost  from  the  incident  wave  via  interference  with  the 
scattered  wave,  and  ~^|gr,  |2  shows  it  is  balanced  by  radi¬ 
ation  over  all  directions.  In  (41),  the  interference  loss  is  bal¬ 
anced  by  specular  scattering  in  the  forward  and  reflected 
directions,  and  by  fluctuation  scattering  over  all  directions. 
Equation  (41)  is  the  same  form  as  (5:62)  for  the  two-dimen¬ 
sional  case;  using  the  Zemike-Prins6  distribution  function, 
we  showed  before5  that  it  reduced  to  the  required  forms  for 
the  uncorrelated  and  periodic  limits. 

From  (39),  we  construct  the  transmission  and  reflection 
coefficients, 

Tu  =  (1  +  C>„  -  C?22  +  C2|  Wl)/A  R2l  =  2C?2l/D, 

1  —  C(,„  +  *22)  +  C2|  |,|  |.  (42) 

The  corresponding  differential  scattering  cross  section  per 
unit  area, pQrl  =:/?c,  7T[K  ] \GrX  |2  is  determined  by 

\Dl2\Grl  I2  =  yrl  |2|  1  -  C>22| 2  +  \pr2l2\CP2,\2 

+  2  Re  p*  pn ( 1  —  Cp\2  )Cp2l.  (43) 

Using  these  and  (30)  we  write  the  transmitted  flux  as 

I T’, , |2  =  1  —  |R2,|2  —  py ,,  sec  #„> 

IT,, |2  =  1  +  4CRe[W,  -  C|  WIKI  -  Cpf2)]/\D  |2, 
111-®  |2  =  ^llU  —  ^F22|2  +  <r22|C^2l|2 
+  2  Re  Oj|(l  —  Cp22  )Cp2i ; 

^2i=  -  (f*i  +  ?12K/2,  a  =  +aa.  (44) 

For  grazing  incidence,  cos  0o—*O  and  C— ►  oo,  so  that 
(7— >0  and  T—*l.  Except  near  grazing,  for  small  \Cp\,  we 
have 

Tzz  1  +  2Cpx ,  +  2C2(p]t  +^21^12)’  ^~2^2,C, 

| T |2=  1  -  |2C>21 12  +  4C Re  ?u s  1  -  \R  \2  - pvt ,  sec  0O, 
Qr\  ~9rl  =c^[K]\pr]\2,  (45) 

where  we  neglected 

R*(?ii  +  IpiiI2  +  Pi  2^21  +  IF21I2) 

-  2ReO»,,alt  +  p2l02l)  =  ^((T2). 

For  small  ling's  g>Re  g,  or  small  W  and  small  absorp¬ 
tion,  we  iterate  (37)  for  g"  sg'  and  neglect  absorption  in  the 
quadratic  terms,  to  obtain 


+urarjg;^>  =g;,  -^^|g;,|2 
=  -  [aa  +  as)/a0, 

a,=^%r[K]\gr0\2, 

|7’|2sl  —  \2C  Img2,|2  —  p((Ja  +  <j,)  sec  d0,  (46) 

such  that  for  W  s \,a,  approximates  <7, .  The  corresponding 
differential  scattering  cross  section  Q^,  ~  W[K  ]q^,  with 
=  c,  l&o  1 2  as  the  value  for  an  isolated  obstacle,  may  be 
used  to  isolate  W\K  ]  from  measurements  essentially  as  in  x- 
ray  diffraction  by  liquids.  For  the  direct  problem  we  use  the 
Zemike-Prins6  p  and  W\K  ]  for  the  planar  distribution  of 
parallel  cylinders;  for  bounded  obstacles,  we  have  explicit 
results  for  restricted  ranges  and  numerical  results  for/?  based 
on  the  Percus-Yevick  or  other  integral  equation  approxima¬ 
tions.9 

For  kb  small  (small-spaced  particles), 

W[K\  =  2T[0]  +0{k2), 

7r[0)  =  W=\+p^[p(R)-\)dR,  (47) 

where  W  is  proportional  to  the  variance  (fluctuation)  of  the 
number  of  particles  in  a  central  region.  We  obtained  W  in 
closed  form  10,1 1  by  using  statistical  mechanics  theorems  and 
the  scaled  particle  approximation  for  the  equations  of  state13 
for  fluids  of  impenetrable  particles.  For  planar  distributions 
of  correlated  parallel  strips  (2^,)  and  disks  (2^"2),  and  for 
lattice  gas  statistics  (3%),  we  showed 

2r0=i-u?,  7rx  =  (\-wf,  ar2  =  (j-rJ?)3,  (48) 

1  -I-  w 

where  w—pv  is  the  fraction  of  the  xy  plane  covered  by  statis¬ 
tical  mechanics  particles  of  width  or  area  v.  The  upper  bound 
on  w  is  unity  for  2P~,  and  2^0;  for  W2  we  use  w  S  0.84,  as 
measured  for  circular  disks.  The  associated  function 
S=w%y  has  a  maximum  Sn  at  w  —  wn  determined  by 
dwS  =  0,  and  so  does  the  corresponding  differential  scatter¬ 
ing  cross  section  per  unit  area  (ocSn).  If  the  particles  are 
nonabsorbing,  the  sum  of  the  coherent  fluxes 
|r2  +  |R  |2  =  1  —pTTa,  sec  d0  has  a  minimum  at  wn;  if 
they  are  absorbing,  and  cra/as  is  small  then  the  minimum  is 
shifted  to 

u?A  =  wn  +  ajas  \&wSn  |>u>n; 

if  aa  /a,  is  not  small  then  there  is  in  general  no  minimum 
with  variation  of  w,  and  cra  dominates  for  all  w.  See  analo¬ 
gous  discussion  of  random  volume  distributions  for  details. 

For  some  purposes  we  may  compare  monolayer  results 
with  equivalents  for  single-path  transmission  through  a  slab 
region  volume  distribution.  Multiplying  numerator  and  de¬ 
nominator  of  2Cplt  =pa0pll/  2cos  0o  by  the  layer  thickness 
d0  (the  separation  of  the  tangent  planes),  we  define  an  equiva¬ 
lent  complex  index  of  refraction  r\  by 

2C> , ,  =  d0  sec  0o=/k  (77-I  )d  =iA , 

d  =  d0  sec  0O,  (49) 

where  p/d0  is  the  number  of  particles  per  unit  volume,  and  d 
is  the  ray  path  along  the  direction  of  incidence.  Similarly, 
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R  =  2C?2x=iFkd  defines  —  Fas  an  equivalent  Fresnel  in¬ 
terface  coefficient. 

Introducing  (49)  in  (45)  and  (46),  we  exponentiate  iA, 
and  work  to  lowest  order  with 

7~e“(l  +  i|*|2), 

|r|J^e-2Im'1(l  -  |*  |2);u  1  —  2  Im  A  -  |R  |2, 
iA  zzi2C  Im  gn  —  p(oa  +  W a, )/2  cos  60, 

|/I  |22s|2CImg2l|2,  (50) 

where  we  would  use  a,  =  ]|Imgn|2  if  the  scat¬ 

tered  were  not  small  spaced.  To  lowest  order  in  k,  except  for 
pressure  release  scattered,  from  g'r>  =  if  r\  + 
with  r  =  kV /a 0  =  tf(k  2,k 3)  and  independent  of  k,  we 

have 

ReA~Pr-kRe^n/2y0, 

ImA^ip^k  /2y0)  [  Im  Re  P  „  )2  ] , 

\R\2^\Rt^2xPrk/2y0\\  (51) 

For  such  scattered,  |R  |2  =  tf{k 2)  dominates  a1  =  tf[k3,k4) 
in  the  scattering  losses  of  a  transmitted  ray,  except  for  special 
angles  corresponding  to  small  & 2\- 

The  limitations  of  (50)  and  (45)  are  indicated  by  the 
complete  expression  T of  (42).  For  detailed  computations,  we 


r-in^=i±f.  |r|>- ,  +  1±J^. 


tan  0T  —  ■ 


2  Im  Z 


Re  T  1  -  |Z  |2 
Z  =  C$ it,,  +  C2^|j^2i/(1  —  Gy22) 

=  c>„  +  C  W2l(l  +  C>22)  +  <?{k4).  (52) 

For  oa  at  least  as  small  as  a, ,  the  transmitted  phase  to  fi(k  3) 
in  terms  of  r,y  =  C  Im  g"  with  g"  sg'  for  small  a/b  is  given 
by 

tan6>r~  2<r|1  ~^r2,T22)  +  -  rl2r2lr22). 


®r~2(r,|  —  t12t2,t22  —  t  /3). 

This  provides  a  more  complete  result  than  in  (51).  Similarly 
from  the  complete  form  of  R  of  (42),  we  specify  the  reflected 
phase  by  Im  R  /Re  R  =  tan  0R  such  that 

jnn  Q _ r2l(l  ~  rnr22  +  rl2T2l) 

r2l(T22  +  T|i)  +  CRe^2| 

_ 1  ~  rllr22  +  r2lrl2 

T22+^l 

(54) 

=(r.i  +  r22Hl  -  r12r2l)  -  (r?,  +  ri2)/3  +  tt/2, 

where  we  dropped  Re  ?2,  for  negligible  losses.  The  set  of 
four  phases  corresponding  to  T,  R,  T\  and  R '  of  (30)  for 
negligible  5 satisfy  &R  +  -  &T  -  @T.  =  n,  where  &T  • 

follows  from  0T  of  (53)  by  interchanging  1  and  2,  and  where 
=  &n  of  (54).  Thus,  to  the  present  degree  of  approxima¬ 
tion.  the  sum  of  the  reflected  phases  diffed  by  180*  from  the 


sum  of  the  transmitted  phases  as  discussed  earlier  for  the 
lossless  periodic  case  (12:15).  Similarly  for  negligible  losses 
we  require  | r |  =  |7”|  and  |J?  |  =  |R  '|;  this  will  be  demon¬ 
strated  in  terms  of  explicit  results  for  ?  we  obtain  subse¬ 
quently. 

The  general  forms  (39)  and  (42)  can  be  simplified  consid¬ 
erably  for  normal  incidence  k  =  2  and  arbitrary  scatterers, 
or  for  scatterers  symmetrical  to  reflection  in  the  plane  of 
centers  and  arbitrary  k.  For  the  first  case,  with 
k,  =  —  k2  =  z,  the  reciprocity  theorem  (6)  gives  g,,  =  g22, 
and  although  g21  need  not  equal  g12,  we  may  factor  the  de¬ 
nominator  to  obtain 
T=Tn  =  r22  =  i[&+  +  ®_), 


—  =  'R„  =*  =  n^+-^_], 


1  +zd 
1  -Z, 


z±  =C(p±/),  ?  =  ?it=?22,  ?’  =  (?, tfu)1 


\ftl2/ 


„  ?r\  +  ,  ?r\  -  '■?*. 

1551 

This  is  restricted  to  normal  incidence,  but  the  scatterers’ 
shape  and  orientation  are  arbitrary.  If  the  scatterers  have 
inversion  symmetry,  then  gc2x  =  ?l2  —  p  and  a  =  1  for  arbi¬ 
trary  orientation  with  respect  to  z  =  0. 

The  same  simple  forms  as  for  a  =  1  also  apply  for  arbi¬ 
trary  k  for  scatterers  having  reflection  symmetry  in  2  =  0, 
with  no  requirement  of  inversion  symmetry.  For  pr,  = 

+F*),  F$  =/$/(l±Z±), 

A  ~  ?iO  ±  Z±  ~  C/ob  =  C(^oo  ±  Poo'Yj 

T=l+  ~  Z+Z- 

1  -  Z+  1  -  Z_  D 

R  z+  z_  z+-z_ 

1-Z+  1-Z_  D 

D—\—  Z+  —  Z_  +  Z+Z_.  (56) 

The  functions  and  /*"  are  twice  the  components  of  ? 
symmetrical  and  antisymmetrical  with  respect  to  reflection 
of  k  (or  f )  in  the  plane  z  =  0.  Equivalently,  or  f~  are  the 
scattering  amplitudes  for  protuberances  on  a  base  plane  for 
which  dxi/>  or  ^equals  zero,  a  rigid  ( -f  )orfree(  —  (base.  For 
<t> '  incident  on  such  structures  the  corresponding  coherent 
reflected  waves  are  ±  ^(1  +  Z±  )/(l  —  Z±  )  with  Z±  as 
normalized  impedances.  We  discuss  such  problems  in  detail 
in  a  sequel. 


V.  MULTIPOLE  EXPANSIONS 

For  various  detailed  applications  or  computations,  we 
reduce  G  [g]  of  ( 1 6)  by  expanding  the  scattering  amplitudes  G 
and  g  in  Fourier  series.  This  provides  algebraic  systems  for 
the  corresponding  scattering  coefficients  A  and  a  in  terms  of 
distribution  integrals  (continuum  analogs  of  the  lattice 
sums  for  the  periodic  cases).5,6  In  Ref.  7,  using  the  general 
form  G  [g]  we  derived  the  general  algebraic  systems  A  [a]  for 
arbitrary  configurations,  and  indicated  how  the  same  results 
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followed  from  more  detailed  manipulations  starting  with  ex¬ 
pansions  of  all  fields  in  terms  of  special  functions.  We  includ¬ 
ed  essentials  of  special  function  derivations  for  the  planar 
periodic  cases3,4  to  stress  several  points,  but  now  we  start 
directly  from  (16). 

A.  Cylinders 

For  cylinders,  we  use  the  trigonometric  series 

grO  =  X  a"e‘ne’ 

ft  =  —  oo 

0)=  f 

m  =  —  oo 

<?*>  =  X^e"10, 

where  the  angle  independent  coefficients  satisfy 
=  (—  1 )" +  ma  _  m  „  (from  the  reciprocity  theorem).  Substi¬ 
tuting  into  ( 16)  we  obtain  in  general 

An  =  a„  +  -  v  (57) 

vm 

and  for  circular  cylinders  as  in  (5:93) 

An=a„[e-^+lAm^m_v],  (58) 

where  the  an  are  independent  of  angles.  The  distribution 
integrals  3Pn  =  (  —  1  )"«#'’_  „  are  given  by  (5:93)  ff, 

=S{e-'v,e‘  +e-,",’r-s*1) 

=  dxp(x)Hn(kx)W^B-[  -  1)" 


with  H„  —  H™  —  J„+  iN„  and  corresponding  decomposi¬ 
tions  =  f  +  or  S  =  Re  S  +  Sc.  The  second  form  in 
(59),  which  follows  from  the  first  by  applying  the  spectral 
integral  in  S  to  the  special  cases  of 

=  rHn(kr)ein 6 

at  hand,  makes  the  dependence  on  k  and  0O  explicit.  [In  Ref. 
5  (1953),  we  derive  the  second  form,  as  well  as  (1 1),  by  ex¬ 
panding  all  fields  in  terms  of  special  functions.] 

For  even  and  odd  n,  we  decompose  (59)  into  two  sets 

S^2n  =  2S  cos  2nOc 

=  2  I  dx  p(x]H2n  (kx)  cos  (kx  sin  0O), 

J  (60) 
jr2„  + ,  =  —  2/S  sin  (2«  +  1)0C 

=  -2 ifdx  p(x\H2„  + ,  (kx)  sin  (kx  sin  0O). 


\“/b  |  ir  _  -  1 

D  =  b-b  =  [\-w)/p,  (62) 

where  [u]  is  the  closest  integer  to  u  from  below.  We  may  also 
use  the  residue-series  representation3 

PP(x)  =  X  ’  1  +Dyv=e-r’b; 

v--«  b  +  De **■ 

yv  =  r±  H  =  -  «w  ±  tf,vi ; 

derived  by  taking5  the  Laplace  transform  of  (62),  and  then 
the  inverse;  see  Ref.  5  for  simple  approximations  of  the  roots 
near  the  sparse  (u>=;0)  and  periodic  (to  ~  1 )  limits.  From  geo¬ 
metrical  considerations,  or  directly  from  (62),  we  obtain  the 
virial  expansion  in  powers  of  w  =  pb,  e.g.,  to  t?(w2)  in  terms 
of  u  =  x/b, 

p  =  0  for  u  <  1,  p  =  1  +  w{2  —  u)  for  l<u<2, 
p  =  l  for  n  >  2.  (63) 

The  Fourier  transform  ofpp(x),  except  near  the  periodic 
limit,  is  given  by 

P  [Kc  ]  =  [Kc  ]  +  (2irp/k  (6(sin  8C  -  sin  0O); 

3>{JT)  =  1  +  ^(JT) 

=  l+2inf  [p  —  1) cos-^Tn du  =  (1  + 

164) 

3)  _  w  sin  JT  nr2)!  —  cos  JT) 

2  ~(l-w)JT+  (1  -w)2JT2  ’ 

Jf  —  kb  (sin  0  —  sin  80), 

where  —  3>  is  the  transform  of  the  direct  correlation  func¬ 
tion.  See  Zemike  and  Prins6  for  the  original  derivation  of  W 
and  for  plots  of  W  and  pp[x)  for  several  values  of  b/b  —  w. 
For  small  kb,  or  for  0tz0 o,  from  (64)  by  expanding  sin 
and  cos  or  from  the  virial  expansion  (63), 

W\3T)  =  W  +  BA2  +  i?(  JT4),  W  =  [\  -  in)2, 

B=  5rw(l  -3in/4)/3,  (65) 

with  W  =  as  in  (48).  For  large  JT(1  —  tv),  from  (64)  in 
terms  of  3> ,  or  from  integration  by  parts  and  using  the  scaled 
particle13  p(b), 

5T(^T )  ~  1  -  [  2t v  sin  JT/(  1  -  w)JT  ]  (66) 

and  W  tends  to  unity  in  an  oscillatory  fashion  as  in¬ 
creases.  We  also  showed5  that  in  the  periodic  limit 
b—*b,w—+ 1, 


For  normal  incidence  (0O  =  0),  the  only  set  that  arises  is  _  i  +  ( lUL\  y  $(Sin  $c  _  sjn  $n ) 

/-  V  kb  )  h  «,  —  oo 

^ta=2j  dxpWM-  (61)  sin  0n  =  sin  0O  +  2  nir/kb, 


See  earlier  work3  for  discussion  and  reduction  of  the  same 
algebraic  systems  (57)  and  (58)  in  terms  of  the  different 
(lattice  sums)  appropriate  for  the  periodic  case. 

The  2T%  are  determined  by  the  Zemike-Prins  pair 
function5'6 


sin  0n  =  sin  0O  +  2  nn/kb, 

corresponding  to  the  grating  of  spacing  b.  The  present  paper 
considers  only  the  low-frequency  periodic  case  kb (l 
±  sin  0O)  <  for,  the  case  of  one  propagating  mode. 

We  obtain  f  n  from  the  first  form  in  (59)  by  replacing  S 
by  Re  S, 
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cf  n  —  fn  +  f'n  ~  $nO  > 

fl  =C{e-lne“  +  e-*nirr-a°)}, 

/'„  =\^nr[K]\e-ine  +  e-ln{’r-e>\,  (67) 

where  C  =  p/k  cos  0O  —  {kb  cos  0o)~l,  and  [  J 
=  2  cos  2nd,  —  H  sin(2/i  +  1)0 .  The  specular  contribution 
/^n  is  the  same  as  for  the  one-propagating  mode  periodic 
case3  with  b  replaced  by  b;  the  term  S„ 0  converts  coefficients 
a  ofg  as  in  (62)  to  corresponding  coefficients  a'  of g'.  We  may 
also  write 

fn \=p\dx[p-\)Jn[  } 

with  [  ]  as  in  the  second  form  in  (59).  For  small  kb,  from  (64) 
and  (67),  or  from  f(J)  with  J„(x)x(x/2)n  /nl,  we  have 

f'o  ~  IT  +  (kb  )2B  ( 1  +  2  sin2  90)/ 2, 

f\  ~i(kb  )2  Bsin0o, 

f\  a  (kb  )2B  /4;  f’2n  ,/\n  =0(*2").  (68) 

The  next  terms  are  two  orders  higher  in  k,  and  the  result 
fn— f°r  kb—* 0  corresponds  to  small-spaced  scat- 
terers. 

We  construct. Y by  using  —  iSe  in  the  first  form  of  (59) 
or  Nn  in  the  second.  For  the  second  form,  S%pNn  [  ] ,  for 
n  —  0  and  1  we  use  So  (p  —  1 W  [  ]  because  the  correspond¬ 
ing  integrals  without  p  vanish;  for  n>l,  we  use 
Sb  +  SbiP  —  !)•  Correct  to  O  (w2)  and  to  leading  terms  in  k 
for  small  kb 

,ro~mo=  -(2/t r)/0, 

/°  =  2  u>  P(P  ~  l)ln  -2-  </«=/„  In  A-  -  /* 

Jo  ckbu  ckb 

cs  1.781; 

(69) 

/„  =  2wj{p  —  1  ]du  =  —  2iu  +  w2, 

I,  =  2 wj(P  —  1)  In  u  du  s;2u;[  1  4-  u;(21n  2  —  J)]; 


— (2n)(  1  —  — — ~  ^ 

2  w  \  w  J 

+  n(2n  +  \)(--~?)2[£(2n)+£(2n  +  1)J. 

w 2 

If  w— »1,  then  the  asymptotic  result  reduces  tof  (2 n),  the  Rie- 
mann  zeta  function,  as  before  for  the  periodic  case.3  For  the 
dipole,  1 2/ 2  (2)~  ?r76  a  1 .645. 

B.  Bounded  obstacles 

For  bounded  scatterers,  we  use  series  of  spherical  har¬ 
monics 

&o  =  £  ±  <Y:(r), 

n  s=  0  m  —  —  rt 

°:  =  <(io)  =  'Z<i‘Yv-,1(io); 

vn 

G *  =  2A  rm  =  (cos  0 ]e**, 

nm 

p-  "  =  P”(  —  ir  (n  -  m)\/(n  +  m)l  =  PfC?, 
a^  =  (-l)n  +  va-nr-m. 

From  (16),  in  terms  of  coefficients14  d f  of  the  expansion 
Y->Y'=^d, 
we  have  in  general 

A:=^  +  iia^A'syd/'M\t)jn-/t  <7i> 

v/i  st  “  \  V  I S/ 

and  for  spheres 

A :  =  ( -  ir  a.  [  y-  m(f0)  +  2A  ~nm  |')*v m  • 

(72) 

These  systems  are  the  same  as  considered  before4  (p.  644  ff) 
for  the  doubly  periodic  planar  lattice,  but  the  present 
represents  distribution  integrals.  We  have 


-  /(2/ir)/0sin  0O=  -  /(2/tt)u;(2  -  w)  sin  6>0.  =  S[  Y”(ic)  4-  F^)] 


Ifu;~  1,  then/^~ln  27r,and/° - In  (Arr/ckb ).  The  correc¬ 
tions  to  are  0(k2).  For  n>  2,  from 

V  \  X  / 

we  obtain  the  dominant  terms  for  small  kb 

r  -  ^  _  (2/i  —  1)!  (  2  \2\ 

2"~(*6)2"  tt  l kb)  l2"’ 

r.  (70) 

+  i  ~  ~  sin  0O,  I2n  =  2w  -£-du. 

Ji  u2n 

For  small  w,  to  O  (ur2), 

I2  =  2 w[  1  +  w(l  —  In  2)]  =2u^l  -|-  u/).307), 

7„=-i!S-(l+-2-(2»-3  +  2-"-)); 

2n  —  1 V  2n  -  2  7 

the  next  term  in  I2  equals  2u»3(l  —  6  In  2  —  3  In  3)=:2u/6.45. 
For  w  near  1, 


=  (  -  1  )mP  ?(0)pj d  R  p(R)eik'Rinhn  (kR  )elm*  (73) 

with  h„  =  h ««  =jn  4-  in„  and  the  corresponding  decompo¬ 
sition  ^  =  /4- LV  as  before.  Either  form  requires  that 
n  —  m  be  even:  for  odd  parity, 

y:(?)4-r:(n=ir:(f)[i  +  (-i)"-m] 

vanishes,  and  so  does  P  ”(0).  The  second  form  follows  from 
the  first  by  using 


^Y:(ic)^rY:(t]hn(kr) 


and  specializing  to  r\\v,ir  4-  <P )  to  correspond  to  the  dis¬ 
played  e*'R.  For  negative  m,  we  have 

*Y"«  c:n-0). 

No  explicit  analog  of p  of  (62)  exists  for  two-dimension¬ 
al  distributions,  but  for  circular  symmetry  p(R  )  can  be  com¬ 
puted  numerically  from  the  Percus-Yevick  integral  equa- 
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tion  or  from  other  approximations  (hypemetted  chain  or 
Bom-Green)  or  from  machine  simulations.9  The  virial  ex¬ 
pansion  for  identical  impenetrable  circular  disks  of  diameter 
b  analogous  to  (63)  can  be  obtained  solely  from  geometrical 
considerations;  in  terms  of  u=R/b  and  w  =  irb  2/4,  we 
again  have  p  =  0  for  u  <  1,  and  p  =  1  for  u  >  2,  but  now  the 
first  shell  is  specified  by 

(74) 

The  structure  factor  for  circular  symmetry 


2T(JT)  =  1  +  Sw  f  (p-  \\J0(JTu)u  du, 
Jo 

P r/kbf  =  (a-a0)2  +  (p-Po? 


(75) 


=  sin2#  +  sin2  90  —  2  sin  6  sin#0  cos(^>  —  <p0) 


may  be  evaluated  numerically  by  means  of  integral-equation 
approximations9  for  p. 

Using  (74)  and  (75),  we  obtain  for  small  3Y~ , 

2T(  JT)  =1-4  w  +  -^-u/2 
n 


+  <?(w3)  +  J'(k 3) 

s  1  —  4u>  +  6.6159u/2  +  \3ir2w(\  —  5.9746w) 

+  (76) 

The  first  three  terms  of  7P~(Q)  =  Tf'  in  (76)  differ  from  the 
corresponding  terms  of  W2  in  (48)  in  that  W2  approximates 
the  third  by  Iw1.  A  more  complete  expansion  of  W~(3i/') 
would  show  the  same  structure  as  (65)  with  2^(0)  =  W  as  a 
factor.  We  use  W ~  JP~2  and  2?  ~  uj(  1  —  6w)/2  in  the  follow¬ 
ing. 

For  large  Jf ,  we  integrate  (75)  by  parts  and  use  the 
scaled-particle  approximation  for  p[b),  and  Jl{x)~{2/ 
rrx),n  sin(jt  —  zr/4),  to  obtain 

,  4w(2  -  w)(2/n)xn  sin(JT  -  jt/4) 

- (T3^ - 3rrr— ’  (77) 

which  tends  to  unity  more  rapidly  with  increasing  JT  than 
the  corresponding  result  in  (66). 

For  aligned  elliptic  disks  with  principal  diameters  b2 
and  by,  we  write 

R  =  u(xb2  cos  t  +  yb3  sin  r)=R(u) 

with  u  =  1  as  the  exclusion  ellipse,  and  use  d  R  =  b2b3u 
du  dr  to  integrate  over  similar  ellipses.  For  such  elliptically 
symmetric  statistics,  we  work  with  (74)— (77)  in  terms  of 

{JT  /k  )2  —  b\(a  —  ar0)2  +  b\(fi  —  f30f 

and  w  —  ptrbyby/^  on  the  basis  of  geometrical  consider¬ 
ations.  The  geometrical  distribution  of  hard  circular  disks  of 
diameter  b,  and  variable  R  =  «6R  with  u  =  1  as  the  exclu¬ 
sion  circle,  determines  the  valuep(J?  J^=  p[  u  j  corresponding 
to  the  equiprobability  circle  R  =  ubR  with  origin  of  u  at  the 
center  of  one  disk  regarded  as  fixed.  A  homogeneous  defor¬ 


mation  of  space  that  leaves  the  origin  fixed,  a  centro- affine 
transformation  that  converts  the  circular  disks  (b )  to  aligned 
elliptic  disks  (fi^fi^such  that  b2b3  =  b  2),  converts  the  equi¬ 
probability  circle  R  =  ub  R  around  the  origin  to  the  corre¬ 
sponding  ellipse  R(u)  on  which  we  usepf  u ) . 

For  circular  symmetry,  the  second  form  in  (73)  reduces 
to  two  sets  corresponding  to  n,m  both  even  or  both  odd: 


=  P7(0)in  ~  melm*°2irp 

x[piR  )hn(kR  )Jm[kR  sin  0O)R  dR, 

P\™(0)  =  (  -  1)"  -  m(2n  +  2m  -  l)!!/2n "  m(n  -  m)\, 
P^,‘(0)  =  (2n  +  2m  +  1)F^(0), 
n!!  =  n(n-  2)!!,  1!!  =  0!!=1.  (78) 

For  normal  incidence  (60  =  0),  only 

=  P°2J0)(  -  l)"2upj“°p(R  )hln[kR  )R  dR  (79) 

is  required.  More  generally,  for  elliptically  symmetric  statis¬ 
tics,  the  symmetry  properties  of  (73)  are  roughly  similar  to 
those  of  the  lattice  sums  for  the  rectangular  unit  cell  dis¬ 
cussed  before  for  (4:18),  p.,  646.  Thus  there  are  four  sets  of 
J^s,  two  for  even  indices  and  two  for  odd: 

-  (  —  1 )" P !:(0)J dR ph2„ (kR  ([cosA'cos  Tcos  2m0 

—  /'  sin  X  sin  Y  sin  2 m<P  ], 

^2?:ti={-i)nPi::,'\o)jdii 

X  ph2n  +  ,  [sin  X  cos  Y cos(2m  +  1  )0 
+  i  cos  X  sin  Y  sin  (2m  +  1)0  ] ;  (80) 

X  =  ka<yR  cos  0  =  kajo2u  cos  r, 

Y  =  kfiji  sin  0  =  k0gb3u  sin  r, 
where 

J/*oo  f2ir 

d  R  =  6263J  du  uj  dr, 
and 

J»2  it  fir/  2 

dr  =  4  dr. 

0  Jo 

Subsequently,  we  use  2X*  and  to  indicate  components 
involving  cos  m0  and  sin  m0.  For  normal  incidence 
(a0  =  0o  =  O,X  =  Y=  0),  the  only  set  that  arises  is 

jnr=(-irFf:(o)JrfR 

X  phlm(kR  (cos  2m0  =  2^7 ■  (81) 

In  (80)  and  (81)  we  use 

R  =  u(b  \  cos2  r  +  b 2  sin2r)l/2 
and 

el<p  =  u(b2  cos  r  +  ib3  sin  t)/R. 

Essentially  as  for  (67),  for  all  cases  with  n  —  m  even,  we 

have 
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f*  =  2CY  ™(f0)  +  -  8«> 

fr  =  ^W[K)Y:(t).  (82) 

The  specular  contribution  in  C  —  rp/k 2  cos  60  is  the  same  as 
the  result  for  the  doubly  periodic  array4  with  p  ~ 1  replacing 
the  area  of  unit  cell.  For  small  kb  and  circular  statistics,  from 
(82)  in  terms  of  7E\K\  of  (76),  or  from  the  form  f\j)  with 
jn^xn/(2n  +  1)!!, 

f0°^7r  +  (kb)2B{\  +  sin20o), 
f\'^-(kb)2B\  sin^», 

-  (kb  )2B&  f*  ,/‘2\  _  ,  =  0(k 2" ).  (83) 

The  next  terms  are  two  orders  higher  in  kb. 

For  elliptic  statistics,  from  (82)with 

(■W /k)2  =  b\(a~a0f  +  b\(P-P0)2, 
or  from  (74)  and  (80)  in  terms  of  u  and  r, 


-yr°*62~ -311=1^1-  I0~  —  8ui(l  —  £,w), 

B,  =  sf— --L)~0.6413, 

M%  =  (2/nt)K(v),  v  =  1  —  t  ~2; 

- | kb2zz  -  3?)  =/„  sin  e^M00e^  +  M20e~  '^)/2, 
M 2  =  (2/irtv) (2E  —  (2  —  v]K  ];  (86) 

J~n/1 

(1  —  vsin2  T)indr, 
o 

(1  —  vsin2  r)~v 2  dr, 
o 

where  E  and  K  are  the  tabulated  complete  elliptic  integrals, 
and  the  next  terms  are  tf(k ).  The  dominant  contributions  for 
n  >  1  involve 

4,  ~~  8--rt  1  +  wB2  -  wl'2n  (, 

2  n  —  1 


foo-Vr^Bk^blb  +  aD  +  bm+p2)), 
f\  =  -Bk  2](b  2a0  +  ib  \0O), 
f?^-(Bk2/\5)(b\+b\), 
ff-ff-m'Kb2  -62), 
f?  =  t>(k\  (84) 

where  the  corrections  are  tf(k  4).  If  b2  =  b3,  these  results  re¬ 
duce  to  the  corresponding  ones  in  (83).  For  both  sets,  if  k—*0, 
then  for  small-spaced  scatterers. 

We  obtain  by  using  —  iSt  in  the  first  form  in  (73)  or 
nn  in  the  second.  Corresponding  to  (70),  the  dominant  terms 
of  JV  for  small  k  may  be  obtained  by  using 
n„(x)=:  —  (2«  —  1  )!!/jc”  +  1  in  the  form  of  (80).  Proceed¬ 
ing  essentially  as  before  for  the  periodic  analog  (4:90),  p.  653, 
we  obtain  in  terms  of  /  =  b3/b2, 


AS*  2m~_ 

In  ~  , 


(4/,-l)l!(-l)"Pl:(0) 
(kb2)2n  +  l 


B2  =^1- 


/s. — g — r*1 ^ 

tt2  2"-2J,/2  u2"-' 

7;=l<L--iv3),  /;-■ ^V3-iL), 

IT  2  TT  L 

,i=iJcT'~JL):  i  =  te<2  +  v3)=I.M7. 

The  dipole  terms  are  given  by 
^(fc62)73-+5R°2/3  =  ~/2M°, 
jE\(kb2)i/9^\/9  =  I2M\, 

I2  =  8ui(l  +  wB2  —  ui/2)ss8u>(1  +  0.4157m), 

M°2  =  (2/nt  ]E,  M\  =  (2/7r3fv)[(2  -v]E-2(\-v\K]. 

-(87) 


If  t— *1  for  circular  symmetry,  then  Ml  and  M  2— *0,  and  M JJ 
and  M°2—*\\  the  leading  terms  of is  then  tf(k  _l).  The 
dominant  terms  for  circular  symmetry  are  the  subsets 


1 

'/r  In  +  1  * 


ni::,' 


(4 n  +  !)!!(  -  \)"Pl::;(0) 
(kb2)2n+l 


+Mi:+2e-  n 


I0  =  8m J  (p-\)du,  I2n  =  8m du, 

)  =  1  C  Jr 

rr  Jo  (cos2  r  +  t2  sin2r)  +m  + 1/2 


where  the  M ’s  may  be  expressed  in  terms  of  elliptic  integrals. 
The  present  IM  plays  the  same  role  as  the  periodic4  L\  as 
before  we  need  consider  only  t>  1  explicitly,  and  obtain  re¬ 
sults  for  t  <  1  by  using 

Mtfr -')=[- vr^+'M  2:(/). 

For  n  =  0  and  1,  in  terms  of p  of  (75), 


jy*  .  (2 n  -  l)H(4n  -  1)!! 

2n~  2"n\(kb  )2" + 1  2"’ 

,  (2m  +  l)H(4n  +  1)!!  r 

%/r  In +\ 


'n\(kb  )2 


',  (88) 


as  obtained  directly  from  the  form  (78). 

Our  consideration  of  W(ff)  and  for  elliptic  sym¬ 
metry  based  on  p(R)  —  p{u],  with  p\u\  as  the  values  of  the 
radial  distribution  function  for  circular  symmetry,  have  led 
to  approximations  in  terms  of  elliptic  integrals  which  make 
explicit  the  departures  from  maximal  symmetry.  By  inte¬ 
grating  over  symmetrical  ellipses,  we  isolated  integrals  of 
p\u\ over  u  that  are  identical  with  results  for  circles,  so  that 
existing  approximations  and  computations9  for  circles  can 
be  applied  directly. 

In  distinction  to  the  one-dimensional  distribution,  for 
which  the  periodic  case  is  a  proper  limit  of  p,  there  are  no 
comparable  representations  or  approximations  for  a  two-di¬ 
mensional  distribution  function  />(R).  There  is  no  unique  pe- 
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riodic  limit  in  two  dimensions,  and  a  great  deal  more  struc¬ 
ture  (noncentral  potentials)  would  have  to  be  introduced, 
and  more  than  geometrical  arguments  would  be  required,  to 
trace  the  evolution  of  a  gas  of  rigid  disks  through  the  liquid 
to  some  specific  periodic  limit.  For  distributions  governed 
solely  by  the  geometry  of  the  exclusion  region,  a  realistic 
bound  corresponds  to  the  two-dimensional  analog  of  an 
amorphous  solid;  instead  of  using  the  close-packed  value 
w ,  2:0.785  for  a  square  array  or  wh  2:0.907  for  a  hexagonal 
array,  we  use  the  experimental  value  wd  2:0.84  for  the  den¬ 
sest  random  packing  (which  approximates  the  mean  of  ws 
and  wh ).  Nevertheless,  the  earlier  results  for  the  rectangular 
lattice4  with  the  square  as  a  special  case  provide  rough  (less 
symmetrical)  analogs  of  the  present  forms  for  the  ellipse  and 
circle,  as  well  as  estimates  for  bounds  for  the  present  case. 

In  the  explicit  low-frequency  approximations  of 
for  the  rectangular  array  given  by  (4:14),  p.  659  ff,  the  func¬ 
tion  L  (t )  corresponds  to  the  present  IM[t ).  For  the  square 
lattice, 

L°0(\\~  -  3.900,  L°(l)2;9.03,  £|(l)  =  0; 
and  for  rectangles 

L  2(2)2:  —  2.521,  L°(2)^4.049,  1^(2)2:1.856; 
Z2(4)2:  -  1.135,  L“(4)^2.815,  Z*(4):=2.267. 
Normalizing  L  ™{t )  by  dividing  by  Z?(l),  we  write 
m2  (2) -0.646,  m§(2)  =  0.448,  m^(2)  =  0.206; 
m°  (4)2:0.291,  m§(4)  =  0.312,  m^(4)2:0.251. 

The  corresponding  values  for  the  present  problem  are 
A/°„(  l)*U#5(l)  =  0,and 

Af  2(2)2:0.686,  4/2(2)^0.385,  4/  2(2)  2: 0.06 16; 
4/2(4)2:0.446,  4/2(4)2:0.171,  4/^(4)2:0.0447.  (89) 

The  normalized  monopole  terms  are  larger  for  the  ellipse 
than  for  the  rectangle,  and  the  normalized  dipole  terms  are 
smaller,  i.e.,  M0(t)>m0{t)  and  4/2(r )  <  m2(r ).  The  ratios 
4/  2  /A/  2  are  smaller  than  m\ /m° ;  with  increasing  t,  the  first 
tends  to  ( and  the  second  to  1. 

VI.  CLOSED-FORM  APPROXIMATIONS 

The  algebraic  systems  (57)  and  (71)  may  be  expanded  in 
series  form,  truncated  and  solved  in  closed  form,  or  investi¬ 
gated  numerically.  We  illustrate  truncation  by  retaining 
only  the  monopole  and  dipole  terms  for  elliptical  scatterers 
and  statistics.  The  results  are  applied  to  obtain  low-frequen¬ 
cy  approximations  for  small-spaced  scatterers  with  empha¬ 
sis  on  dominant  multipole  coupling  effects. 

The  isolated  scattering  amplitude  has  the  form 

g«  =a0  +  M-k  =  a0+  £ 

I 

=«o+  i  x  wrf' 

t-u-i 

r,  =  cos  9,  Yi  —  sin  9  cos  <p,  y3  =  sin  9  sin  <p\  (90) 

where jk  is  a  symmetrical  dyadic  with  principal  values  a, 
along  and  n  =  3  or  2  for  ellipsoids  or  elliptic  cylinders. 
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For  the  corresponding  multiple  scattering  amplitude,  we 
write 

G*  =A0  +  r-A  =  A0+  X^y,- =  X^.y,.  y0=T 
1=1  1=0 

(91) 

We  work  explicitly  with  n  =  3,  and  set  y3  and  <p  =  0  for 

n  =  2. 

From  the  algebraic  systems  (57)  and  (71),  or  from  the 
functional  equation  (16),  we  obtain 

A0  =  a0(  1  +  +  ^-A),  A  =  a-(k  +  3t. A0  +  «5p.A); 

02  +  y3? 03,  |92) 

&  =  3TX  ,zz  +  JT22xx  +  JTnxy  +  ^32yx  +  ^33yy, 

where  the  JP5,,  =  ^ ,  are  linear  combinations  of  the  distri¬ 
bution  integrals  we  considered.  For  the  planar  array  of  cylin¬ 
ders,  in  terms  of  (59)  ff, 

^oo  =  ^0.  ^n  =  \(^0  +  sr& 

<^22  =  J  (^0-^2).  (93) 

For  bounded  obstacles,  in  terms  of  (73)  ff, 

qu?  9^0  9^1  c  _  oj^ls 

0Z  00  —  ^  0 >  0*  02  —  0*  1  »  0i  03  —  0*  1  * 

xru  +  2^2), 

QC/?  _  1  7s  _  1  /  36*50  36*50  \  1  1  Q '/?lc 

23  —  2  »  22  ““  i(^0  ^2  I  <  i**  2  » 

<^33  =  (94) 

The  superscripts  c  and  s  correspond  to  terms  in  cos  m<P  and 
sin  m0  as  indicated  after  (80).  Thus,  in  general,  the  four  A ’s 
are  specified  by  four  simultaneous  equations  in  terms  of  the 
a’s  and  !%*%.  First  we  consider  the  simpler  systems  that  arise 
for  several  special  cases  in  order  to  display  the  structure  of 
the  multipole  coupling;  then  we  obtain  closed  low-frequency 
forms  of  (92). 

If  the  principal  axes  of  the  dipole  are  along  xj/j  (the 
symmetrical  monolayer),  then  (92)  reduces  to  two  sets. 

At=a,{tf  +<r„ Aty,  A,  =o,(y?  +  X*V/)J 

i*  1 

i  or  j  =  0,2,3,  (95) 

corresponding,  respectively,  to  the  decomposition 
G  =  Ga  +Gs  into  components  antisymmetrical  and  sym¬ 
metrical  with  respect  to  the  plane  of  centers  ( z  =  0).  Intro¬ 
ducing  the  self-coupling  coefficients, 

we  write  the  antisymmetrical  component  as 
A\  =a,y?/(l  —  n)  =  j/  ,y?, 

G^(?.k)  =  ^,y,  =  j/.y.rf,  (96) 

and  specify  the  symmetrical  component  Gs  —  A0  +  A2y2 

+  A  iYj  by 

A0  =  s/  o(l  +  A\3f> 02  +  A-sJfi? 03), 

A2  —  s/  2  (y?  +  A0 20  +  Ay^if 23), 

Aj  —  sly?  +  Aoft* 30  +  A-i&f 32).  (97) 

For  normal  incidence,  y?  =  y?  =  0,  and  the  displayed  lat¬ 
tice,  sums  vanish;  they  are  not  in  the  subset  (81).  Then 
Gs  =  sf  0  and  G  =  s/0  +  s/ , Ylyfii- 
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For  incidence  in  the  zx  plane,  only  the  arise  and 
only  3fa2  =  2f\c  survives  in  (97).  Thus  Gs  =A0  +  A2y2  in 
terms  of 

A0  =  £/0(l  +  A23f 02),  A2  =  s/  2(7^2  +  20); 

A0  =  rf0{l  +J^02s/2f2)/D, 

A2  =  s/2(f2  +  jr2o-*/o)/A  D  =  1  -  s/0s/2{^20)2, 

and  the  multiple  scattered  amplitude  equals 
G(r,k)  =  G„  +GS  = 

-3  Q  +  3 2Y2/2  +  d 0^ 2^02(72  +  7*2  ) 

1  -  rf0s/2(jr02)2 

s/^aAl-a,^,).  (98) 

Here  the  s/0s/2  terms  correspond  to  monopole-dipole  in¬ 
teractions.  This  same  structure  also  arises  for  spheroids  with 
symmetry  axis  along  z  and  arbitrary  k,  and  for  the  two-di¬ 
mensional  problem  (y2  =  sin  9 )  of  elliptic  cylinders. 

Although  the  general  case  in  (97)  is  a  simple  system  that 
can  be  solved  directly,  in  order  to  delineate  cross  coupling 
between  s/2  and  srf  3,  we  first  incorporate  such  effects  in 
coefficients  3) , <6 .  We  have 

A2  —  3  23  +  A0V  23,  A  3  =  ^32  +  Aq^j  32, 

38  22  =  s/  2  (y?  +  y%  s/  3J%f’2iVD, 

<6  23  =  ^2(^20  +  s/3^23^30)/D, 
D=l-^2s/3(Jr23)2,  (99) 

and  38  32,^32  follow  on  interchanging  2  and  3.  If  the  mono¬ 
pole  is  negligible  (e.g.,  if  the  relative  compressibility  of  the 
particles  is  unity),  then  A2  =  38 2}  and  yf3  =  38 32.  From  (99) 
and  (97). 

4  fi—  oil  it  ^ 03^23)  (100) 

\  ^ m'tS 23  +  033>  32) 

and  A2  and  A3  follow  from  (99). 

An  alternative  simplification  of  (92)  is  obtained  by  spe¬ 
cializing  the  forms  for  arbitrary  alignment  to  £3  =  y  and  nor¬ 
mal  incidence  k  =  z.  Only  the  are  nonvanishing,  and 
from 

A0  —  a0(  1  +  Ao&’oo), 

A |  =  fln(l  4-  33* uAt)  +  ai2<^ 22A 2, 

A2  —  02|(1  +  \ \A  |)  +  a 22^2*4 2,  (101) 

we  have 

A0  —  s/  0,  A  |  =  [flu  —  (flna22  —  22]/^’ 

A2  =  a2Y/D\ 

D  =  1  —  Q  |  1 1  —  Q223f  22  +  (®  1 1^22  —  ®?2  j  22! 
G  —  A0  4-  Aty{  +  A2y2.  (102) 

the  same  results  hold  for  the  two-dimensional  problem.  Indi¬ 
cating  orientation  by  g,  =  i  cos  P  —  x  sin  P,  we  rewrite  the 
dipole  terms  as 


A,  =  (a,  cos2  P  4-  a2  sin2 P  — 

A2=  —  (a,  —  a2)sin  PcosP  /D, 

D  =  1  —  a,(<^,  1  cos2  P  4-  ^22  sin2  P ) 

—  a2(^i  1  sin2  P  4-  ^22  cos2  P ) 

+  fllfl2^’n^22-  (103) 

The  forward  scattering  amplitude  for  this  case, 
G2!  =A0  +  A,  with  3fu  =  Jr„(zJ,  is  substantially  different 
from  G  of  (98)  evaluated  for  r  =  k  =  z  cos  P  +  i  sin  P,  i.e., 
from  Gkk  =A0  +  A,  cos2 P  +  A2  sin2 p  with  tj 
=  3? ij  (k).  Although  the  orientation  of  a  single  ellipsoid  with 
respect  to  an  incident  ray  is  the  same  in  both  cases,  the  orien¬ 
tation  of  the  array  is  not,  so  that  the  coupling  processes  dif¬ 
fer.  In  general  (103)  is  considerably  more  complicated  be¬ 
cause  it  is  much  less  symmetrical. 

More  generally,  we  reduce  (92)  to  the  form  G  \g]  of  (32) 
which  we  have  already  analyzed,  and  show  the  relationships 
of  g,g',g",  and  g  explicitly.  Substituting  the  first  form  ofg  of 
(90)  and  the  analog  for  g'  into  (36),  we  obtain 

g'r o=0(i+  fa'*,  «£=(!+  «o)~  V 

a'  =  (I4-  a/n)_l-a.  (104) 

The  amplitude  g'  for  the  radiationless  obstacle,  and  the  coef¬ 
ficients  a'  are  imaginary  for  lossless  scatterers.  In  practice, 
we  start  with  g‘  or  an  approximation  (derived,  e.g.,  from 
perturbation  of  potential  theory  results)  and  use  the  inverse 
relations  a0  =  <*o(l  -  a^)~  \a  =  §'•(/  -  a 7n)~ 1  to  construct 
g.  To  analyze  (92),  we  use  (104)  and  the  decomposition 

/%  =  2Cy?y?, 

(105» 

where  m0  =  1,  and  m,  =  1/n  with  n  =  2  or  3  if  0. 

Substituting  (105)  into  A0  of  (92),  we  make  the  specular 
contribution  explicit 

A0  =  a0[l+(2C-\+^]A0 
+  2C(f1A2  +  f3A3)  +  Jr-A] 
and  then  suppress  —  1  4-  33”oo  by  the  form 

A0  =  { [  1  +  2C40  +  C(k  +  k>A]  +  ^-Aj 

=  A[,l[(l  +  CG00)  +  CG(y0]+<?r-A},  (106) 

jrf’  = _ ^2 _ 

0  1— a0(— 1+^) 

a0  Oq 

1  —  Voicf  60  +  L/^oo)  1  —  alf  w 

Here  a„  =  /(l  —  a,oLAroo)=ui,o30 &  such  that  3>0  is  real 

when  a'0  is  imaginary,  is  the  monopole  coefficient  of  the  ra¬ 
diationless  distribution  amplitude  g'. 

Similarly  for  A  of  (92),  we  incorporate  —  Sy/n  directly, 
and  make  explicit  the  specular  terms, 

A  =  a'-fk  +  C(k  +  k'Mo  +  C(kk  +  k'k>A] 

-H  *'■{<&* A0  +  JP-A) 
and  then  suppress  iP1  to  obtain 
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A  =  ^.j  [k(l  +  CGqq)  +  k'CG00  ]  +2TA0 j, 

=  (I  -a'-Jr'r'-a'  =  (107) 

where  a'  =  (I  —  A'\.  V)~  ‘4'  is  the  dipole  coefficient  of  the 
radiationless  g".  The  brackets  in  (106)  enclose 
((f'>  )  +  (V  <  )  )/2  evaluated  at  z  =  0,  i.e  the  mean  coher¬ 
ent  field  (excess  pressure);  and  similarly  the  brackets  in  ( 107) 
enclose  V((l/'>  )  +  (¥<  ))/2 ik  evaluated  atz  =  0,  the  mean 
gradient  (normalized  velocity). 

We  eliminate  the  cross  terms  in  (106)  and  (107)  corre¬ 
sponding  to  monopole-dipole  coupling,  and  obtain  in  terms 
of  L  =  1  -f-  CG0 o,  M  =  CG0.0 , 

A o  =  (a0  +  a-k)L  (a0  a‘k']M, 

A  =  [a  -f-  a-k]L  -+■  (a  +  a‘k']M, 

a0  =  (i  - 

a  =  rf-(I  -  ajrzr-tf), 

a  —  -s/ a0  =  a0,  (108) 

as  well  as  a  =  (I  —  and  a  — 

Substituting  into  G#  =  A0  +  r-A  yields  G^  = 

+  M  as  in  (32)  with 

=  «0  +  «•(?  +  ?,)  +  ?•«•?,;  r,  =  k,k'.  (109) 

The  three  coefficients  a0,  a,  and  a  hold  for  both  k  and  k' 
because  only  their  tangential  ( xj> )  components  are  involved. 
If  we  replace  k  by  —  k,  then  a0  and  a  are  unaltered,  but 
a(  —  k)  =  —  «(k).  The  equivalent  scatterer  from  which  we 
obtain  G  by  incorporating  the  specular  effects  via  (32),  satis¬ 
fies  the  same  reciprocity  relation  as  G,  i.e., 
^(fyf2)  =  ?[  —  ?2,  —  r ,)  with  either  r ,  or  f2  equaling  either  k 
or  k'.  Substituting  ^  and  ^  of  (109)  into  (39),  (42),  ff, 
provides  explicit  illustrations  of  the  general  considerations. 

The  procedure  (106)-(109)  serves  to  delineate  the  rela¬ 
tions  between  the  various  coefficients  and  amplitudes,  and  to 
display  the  role  of  the  coherent  field  in  the  monopole  A0  and 
of  its  gradient  in  the  dipole  A.  However,  we  derive  ?  more 
directly  from  (32)  by  substituting  <?rX  =  U„  +  f-U  and 
grl  =a0  +  r  a*r  |,  and  denoting  the  corresponding  distribu¬ 
tion  integrals  generated  by  S'01  as  J%^0>  =  df  —  We  ob¬ 
tain  (92)  in  terms  of  \l0,Vl,^0) ,  and  corresponding  analogs  of 
( 106)  and  ( 107)  devoid  of  specular  contributions  (no  terms  in 
C).  This  leads  directly  to  U0  =  s/q(1  +  ^-U)  and 
U  =  j/-(f,  +  Jf^Uo),  and  after  eliminating  the  cross  terms, 
to  U0  =  a0  +  a-x ,  and  U  =  «  +  *  •?,  with  coefficients  as  in 
(108). 

For  small-spaced  scatterers,  and  a'0  and  a'  of  order  kn , 
we  use 

(no) 

where./l7',,,  obtained  from  (70)  and  (85)  via  (93)  and  (94),  is  of 
order  k~”  for  />  0,  and  k  n—*0  as  k—* 0.  Thus 

?r\  ~aO  "+■  aO~Q<)/{l  °0 

«  =  [I  —  a'(3T/«)  -  ii’JT]  ~ '4' 
==(/-i'5r/«)-,4'=i''  +S'4'5r-/n,  (111) 

a'  =  (I  -  ’4'  =  &  ~  ’4'. 

For  lossless  particles,  i'  and  i'  are  imaginary,  i.e., 


&  =  I  —  is  real;  the  factor  &  ~ 1  introduces  a  fc-inde- 
pendent  coupling  correction  to  the  stripped  dipole  coeffi¬ 
cient  a'.The  operator  (I  —  a "W/n)~'  reinstates  radiation 
losses  corresponding  to  fluctuation  scattering  governed  by 
the  packing  factor  W. 

To  illustrate  the  above  we  start  with  Rayleigh’s  re¬ 
sults15  for  elliptic  cylinders  and  ellipsoids  in  the  form 

a^irV,  <€=C-\\ 


a'  =  -  iT/3/(i  +  #?,.)=  -  ins;,  ®=b-  i- 

r  =  K2r-/A,  k'r/Air,  ImC>0,  ImR<0.  (112) 

Here  V  is  the  particle’s  volume,  C  and  B  are  its  relative 
parameters  (compressibility  and  inverse  mass  density  in  the 
simplest  cases),  and  qf  is  the  elliptical  depolarization  inte¬ 
gral16  normalized  to  2g,  =  1;  in  two  dimensions,  in  terms  of 
the  principal  diameters  </,  and  d2  along  z  and  x, 
qy  =  1  —  q2  =  d2/(dx  +  d2).  For  rigid  scatterers,  C  =  B  =  0, 
and  Oq  ~  —  iFt  a\  /(l  —  ?,)•  More  generally  than  in 
(112),  we  replace  B  by  a  tensor  parameter  B  =  (BtJ ). 

If  the  principal  axes  are  along  the  Cartesian  axes,  then 
for  cylindrical  dipoles,  in  terms  of  /2  of  (70), 

i  —  1  +  33'e,,  e,  =  e  =  TI2/2trb 2,  e2  =  —  e. 

(113) 

For  ellipsoids,  in  terms  of  I2M  2  and  I2M  \  of  (87), 
el  =  e,  e2=  -{(e  +  e').  i  (f  ~ 

e  =  ri^'l/Anb e'  =  Sr-I2M\/4nb  \  (114) 

for  both  cases,  =  0,  and  consequently  a"  has  the  same 
form  as  a', 


Qi  =q,+e 


-irsA, 


i  + 

/.  =  1, 


(115) 


such  that  the  packing  effects  are  incorporated  in  new  de¬ 
polarization  factors  Q, .  The  relation  a„  =  =  iPiS  plus 

Eq.  ( 1 1 5)  determine  an  ellipsoid  having  the  same  volume  and 
physical  parameters  as  the  original,  but  a  different  shape. 
The  equivalent  ellipsoid  is  flatter  along  the  array  normal  and 
broader  in  the  plane  of  centers,  the  elongation  being  greater 
along  the  smaller  diameter  of  the  exclusion  ellipse  (along  b2 
for  b2  <  b3).  Corresponding  to  ( 1 1 1 ),  we  have 


-  2rr2(^2  +  is2r,y?/n),  (ii6) 

where  we  keep  only  the  leading  terms  in  T 2,  and  use 
# 2s(Re  ^ )2  | ^ | 2,  etc.  Substituting  (116)  into  (50)  and 
(51)  we  generalize  the  earlier  results  by  replacing  corre¬ 
sponding  to  g’  by  the  present  factor  for  g"  and  thereby  in¬ 
clude  the  dominant  multipole  coupling  effects. 

More  generally,  for  an  ellipsoid  of  arbitrary  orientation 
with  shape^  dyadic  q  =  2^,|,|,  and  dyadic  parameter 
B  =  SRjUiH,  ,  we  have 

a  =  -  tr(i  +  J?.q)-  ‘.^=  -  iris',  &  =  B  - 1. 

(117) 

Writing  the  array  coupling  factors  ef  with  2f,  =  0  as  the 
array  dyadic  c  =  (where  t,  =  £,£,>►),  we  use 
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-  at  Jr  =  »'•€,  and  (I  -  A'Jh~ =  (1 4-  S'-cp  '4'  to 

obtain 

i*=-  /r(I  +  -  ir%,  Q  =  4  +  i; 

«  =  -  iT8  -  ( arr  V«)8-»,  ( 1 1 8) 

with  Q  as  the  net  depolarization  dyadic.  The  corresponding 
equivalent  scattering  amplitude,  obtained  by  substituting 
(118)  in  (111), 

^/rf*?  -  f.».k)  -  arr  2(<^2  +  ?-88k /«),  (i  19) 

as  well  as  the  special  case  (116),  can  be  used  in  (44)  ff  to 
construct  T,  R,  and  \G^  |2  explicitly. 

We  apply  (119)  to  illustrate  statements  for  (52)  ff  with 
reference  to  theorems  for  T  =  Tw,  T'  =  Tva ,  etc.  It  suffices 
to  consider  coincident  principal  axes  for  the  parameter  (B) 
and  shape  (q)  of  the  particle,  and  particle  orientation  speci¬ 
fied  by  g3  =  y,g,  =  z  cos  /?  —  i  sin  /?,£2  =  y  x|,.  For  inci¬ 
dence  and  observation  in  the  zy  plane,  we  require  only 

8  =  (I  +  8'.e)-'.8'=  [(«;,  +s;s;e2)zz 

4  33'12(zi  +  iz)  +  (3S;2  4  33;33;f,)ii]/A 
D  =  1  +  33J  i  Cj  +  S22e2  +  53J  33;e,€2) 

as;  =  a,{i  +  =b,  - 1,  (120) 

where 

as;,  =  as;  cos20+as;  sin2# 
as;2  =  as;  sin2/?  +  as;  cos2/?, 

9S;2  =  SB',  =(3S; -33;)sin0cos#  (121) 

Equivalently,  in  terms  of  Q  =  q  +  e  =  (Q,t)  =  (qtj  +  e,8u) 
with  qi}  expressed  in  terms  of  £,  and  q2  by  the  relations  in 
(121),we  obtain  from  the  form  S3  =  (I  + 

8  =  [(^n  4  ^,^2£?22)zz  +  (^2,  -  ^,^2C21) 

X (zi +  iz)  +  (^22  4  ^,^20,, )**]/£ ', 

D  =  1  +  38 \\Qn  4  ^  22Qn  +  2^,20,2 

+  & j(Q\  1Q22  —  Qnh  (122) 

with  SB  v  in  terms  of  SB ,  as  in  (121).  From  (120)  or  (122)  with 
8  =  (33,;),  we  have 

88  =  (as2,  +  as?2)zz  +  8,2(8„  +  as22) 

X  (ii  +  ii)  +  (822  +  S3  ?2  )ii,  (123) 

for  which  we  use  (93,y  )2  ==  (Re  93(y)2. 

Corresponding  to  the  generalization  of  C?  in  (51),  we 
write  &  =  @1  +  if  with  and  f  in  terms  of  the  present 
results.  Thus,  for  absorption  of  the  order  of  the  scattering 
losses, 


f2‘8'f,  —  SSnj^y,  4  33  niYi&i  +  O2Y1)  +  SS^^,, 

f2.g.8-f,  =  (as2,  +  822)y2y,  +  8,2(S9„  +  ®22) 

Xir-tPi  4a2y,)4(33|2  +33?2)a2a,.  (125) 

Specializing  to  f2,f,  =  k',k  and  k,k',  in  terms  of  k  =  zy  +  ia 
and  k’  =  zy  +  ia,  we  have 

k'-8-k  =  k-8-k'  =  -23,,^  +  8 22cr2, 

k'-88k  =  k-8-8k' 

=  -  («?,  +  4  (33222  4  33?2)a2.  (126) 

Consequently,  from  R2l  =  2C?2X/D  of  (42) 

Poro  _  R<yo  _  R  _  i  (127) 

Poor  Bay  R 

Thus,  although  the  particle  alignment  is  not  symmetrical  to 
the  array,  the  reflection  coefficients  for  arbitrary  imaged  dir¬ 
ections  k  and  k'  are  identical  in  magnitude  and  phase. 

On  the  other  hand  if  f2  and  f,  both  equal  k  or  k',  then 

k-8-k  =  +  as,22ra  4  3322a2, 

k'-8-k'  =  as,,?/2  —  %n2ya  +  9322a2,  (128) 

and  similarly  for  the  8-8  terms.  Thus,  from  (42)  for  Tw  =  T 
and  Taa  =  T\ 

DT- (1  -  C2y\\)  =  -  [DT'  -  (1  -  C2||^||)] 

=  C  (Poo  ~  Pvor ) 

=  iKWl2ya  -  K4 riT  Re  23, 2 
XRe(8„  +  aS22)ra,  (129) 

so  that  T  and  T’  differ  in  general.  However,  for  negligible 
losses,  ^,2sj/Im^,2>  and  \\p\\  -  pooP<r<r  ~  Poa  «  real;  for 
such  cases,  (129)  in  terms  of  real  quantities  A  and  B  reduces 
to 

rp _ A  4 iB  p, _ A  —  iB  17-12 _ it'I2 _ A2  +  B2 

D  ’  D  ’  1  1  _l  1  ~  |D|2  ’ 

(130) 

so  that  the  transmitted  magnitudes  are  equal  for  k  and  k'. 
This  plus  (127)  illustrates  the  general  considerations  of  mag¬ 
nitudes  after  (54),  and  ( 1 26)-(  1 29)  serve  for  the  related  discus¬ 
sion  of  the  phases,  i.e.,  that  the  sum  of  the  reflected  phases 
(0R  +  0R  =  20r  )  and  the  sum  of  the  transmitted  phases 
(0X  4-  0T )  differ  by  ir  for  negligible  losses.  From  (129),  the 
transmitted  phases  differ  by  &T  —  0r.  s=8K  Re  33, to 
first  approximation. 

The  development  of  (1 12)  excludes  pressure  release  ob¬ 
stacles  for  which  a'0  is  not  of  order  k* ,  and  the  case  of  mono¬ 
pole  resonance.  We  include  these  by  an  alternative  form  for 
the  monopole  #>1  (corresponding,  e.g.,  to  underwater 
sound  incident  on  gaseous  cylinders  or  bubbles), 


Im  Q2,  =  KB P2„  ReQ2,  =  A’(/'2,  + 

K  =  prk/Ay0  =  cr, 

0?2l  =  Reft?  -  f2.«.f,),  f2X  =  Imf#  +  f2.8.f,), 

=  (Re  #)2  +  Re(f2-8)-Re  (124) 

The  angle  dependent  terms  for  r,  =  z  cos  6,  +  i  sin  0, 
=  iy,  +  ia,  are  given  by 


irr 

1  -nsy  ’ 


cs-lAL 

n  ’kf’ 


(131) 


where  =  ln[8 /kc{dx  +  </2)],  and  /is  the  electrostatic  ca¬ 
pacity  of  the  ellipsoid.  For  the  pressure  release  case,  we  let 
C—  oo  in  V  =  C-  1  of  (131),  and  R— 00  in  a,'  of  (112)  to 
obtain 
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°o ~-^r  =  ^L>  -M  = 

Jr  2<f0  q, 

(132) 

where  the  monopole  dominates  for  small  k.  The  form  in  ( 1 3 1 ) 
also  suffices  near  the  first  resonance  for  large  finite 
The  resonant  frequencies  correspond  to  F&  JT  =  1,  i.e., 

kl„  s2ir/?''Wo.  (133) 

for  which  the  isolated  monopole  a0  =  a'0/(  1  —  a'0 )  reduces  to 

nA  =  -  1- 

Corresponding  to  (132)  or  (133),  we  neglect  the  dipole 
and  use 

?^a0~a'0/(\  -  a'0Lrw  -a’07n  =  «J/(1  -  *1, 

flo  —  ao/(l  —  do),  (134) 

intermsof^oo  =  ^o  —  3i0of(69)or^4^  =  9}§/W>2of(86). 
The  corresponding  transmissions  and  reflection  coefficients 
satisfy 

T= - \ - =  l  —2C.t  -  =  1  +  R, 

1-2  C>  1  -  2C> 

|r|2  =  i-|^[2+  (135) 

where  the  last  term  equals  —  p y  sec  0O  as  in  (44).  Near  reso¬ 
nance,  we  write 

a"  =  iF€/{  1  -  7V( X  -  oo)]  =  iT<g/(  1  -  k  2/k  \ ), 


k^=7~h>k^>  *  =  (136) 

(1  - x )  ^  b2 

so  that  multipole  coupling  shifts  the  resonances  to  higher 
frequencies.  At  resonance,  ^0A  =  —  \/W  with 
I^oa  I  >  1*0 a  I  -  Consequently 

T  -  *  R  -  ~2C  I,  ^ 

A  2  C+ar’  A  2C  +  5T  2  c+sr’ 

4T2T 

““’-licTW  1,3,1 

and  A  dominates.  With  increasing  w  in  (48),  W  becomes 
small  (zero  in  the  periodic  limits),  and  TA  sO,  /?A  «  —  1, 
and  y ,  =r0;  thus  such  bubble  screens  approximate  perfect 
reflectors. 
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